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ABSTRACT 

The  use  of  shrinkage  type  estimators  for  the  attrition 
rates  in  manpower  flow  models  is  extended  to  include  the 
limited  translation  James-Stein  technique  introduced  by 
Efron  and  Morris.  The  performance  of  these  estimators  is 
compared  with  that  of  several  "natural"  estimators  on  two 
scales:  the  original  scale  of  rates,  and  the  Freeman-Tukey 
transformation  scale  which  was  chosen  in  order  to  give 
shrinkage  estimators  more  efficacy. 

Generally,  the  aggregate  methods  currently  in  use  are 
outperformed  by  MLE ,  transformed  scale  cell  average, 
James-Stein,  and  limited  translation  James-Stein.  The 
results  among  the  four  were  mixed  when  both  global  and  small 
inventory  cell  figures  of  merit  were  compared.  It  is  felt 
that  better  data  aggregation  will  permit  limited  translation 
to  excell  in  low  inventory  estimation. 
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I.  INTRODUCTION 


A.  PURPOSE 

This  is  a  continuation  of  a  pilot  study  started  by  Major 
D.D.  Tucker  in  a  thesis  [Ref.  1]  submitted  at  the  Naval 
Postgraduate  School  in  September  1985.  The  reader  is 
referred  to  Tucker  [Ref.  1]  for  most  of  the  background 
information,  including  a  detailed  discussion  of  the  Marine 
Corps  officer  attrition  and  promotion  structure,  the  officer 
manpower  planning  process,  and  the  attrition  rate  models 
explored  by  Major  Tucker.  Only  information  that  is 
immediately  pertinent  will  be  repeated. 

B .  BACKGROUND 

The  United  States  Marine  Corps  has  about  20,000  offi¬ 
cers.  These  can  be  cross  classified  into  40  military  occupa¬ 
tional  fields  (OF),  31  lengths  of  service  (LOS),  and  10 
grades,  or  12,400  categories  for  manpower  planning  purposes. 
About  half  (6149)  of  these  categories,  called  hereafter 
cells,  are  unoccupied  for  structural  reasons,  e.g.,  due  to 
policy  decisions  concerning  minimum  and  maximum  lengths  of 
service  for  each  grade,  and  the  allowable  grades  in  each  OF. 
Estimates  of  the  attrition  rates  from  these  cells  support  a 
number  of  Marine  Corps  models,  and  accurate  prediction  of 
the  rates  is  basic  to  effective  manpower  utilization. 

The  goal  of  this  pilot  study  is  to  find  efficient  ways 
to  estimate  attrition  rates  (i.e.,the  rate  of  leaving  the 
service,  not  of  changing  OF,  LOS,  or  grade)  for  the  officer 
OF/LOS/grade  categories.  This  is  a  difficult  problem 
because  of  the  large  number  of  cells  with  low  inventory 
figures.  We  will  refer  to  this  as  the  "small  cell"  problem; 
it  is  this  small  cell  problem  that  is  of  greatest  concern  to 
the  builders  and  users  of  these  manpower  models. 


Because  of  the  very  large  number  of  cells  and  their 
heterogeneous  nature,  it  is  wise  to  collect  cells  into  major 
groups,  or  aggregates,  which  conform  to  certain  assumptions 
concerning  their  statistical  behavior.  If  these  assumptions 
are  at  least  approximately  valid  for  the  aggregates,  then 
certain  theoretical  models  can  be  used  to  predict  attrition 
rates  more  accurately  than  current  practices  allow.  Amin 
Elseramegy  [Ref.  2]  explored  the  aggregation  problem  using 
the  CART  routine  with  encouraging  results.  But  these 
results  were  not  available  or  usable  in  a  timely  fashion  to 
be  included  in  the  present  study.  This  thesis  and  Tucker 
[Ref.  1]  assume  valid  aggregations  can  be  found,  and  explore 
the  performance  of  likely  estimation  schemes.  However,  the 
aggregates  used  in  these  pilot  studies  conform  to  current 
Marine  Corps  practice.  These  were  selected  on  grounds  that 
conform  to  organizational  and  operational  considerations, 
and  are  unlikely  to  be  related  to  a  choice  made  on  the  basis 
of  the  statistical  modeling  behavior. 

Current  Marine  Corps  practice  places  all  OF ’ s  in  four 
categories:  aviation  (OF  72,  75),  combat  support  (OF  13, 
25,  35),  ground  combat  (OF  03,  08,  18),  and  other. 
Aggregates  are  formed  from  these  categories  by  taking  data 
by  grade.  Past  attrition  rates,  from  1977  to  the  present, 
weighted  by  subjective  judgement,  are  used  to  predict  future 
attrition.  An  average  attrition  rate  (the  grand  mean)  is 
computed  for  the  entire  aggregate,  and  this  rate  is  used  as 
an  estimate  for  all  cell  attrition  rates  in  the  aggregate. 
We  expect  to  improve  substantially  over  this  method. 

It  should  be  noted  that  the  aviation  category  used  by 
Tucker  [Ref.  1]  included  only  OF  75.  For  continuity,  this  is 
continued  in  the  present  work. 

C.  PROGRESS 

Tucker  [Ref.  1]  showed  the  James-Stein  shrinkage 
estimator  [Refs.  3,4]  can 

(1)  greatly  improve  on  current  methods, 


(2)  improve,  in  a  global  sense,  over  maximum  likelihood 
methods,  and 

(3)  provide  estimates  for  those  small  cells  which  have 
had  no  attrition,  i.e.,  those  cells  whose  MLE  must 
equal  zero. 

Present  work  continues  this  investigation  in  that 

(1)  global  measures  of  efficacy  (risk)  are  decomposed  so 
chat  the  effects  of  the  method  can  be  separately 
examined  for  small  and  large  cells,  and 

(2)  a  class  of  extensions  of  the  James-Stein  estimation, 
called  limited  translation  shrinkage  estimation,  is 
applied  to  the  aggregates  studied  by  Tucker  [Ref.  1] . 

The  main  purpose  is  to  sharpen  the  treatment  given  the  small 

cells . 

D.  RESULTS 

It  appears  the  limited  translation  technique  adds  to  the 
efficacy  of  the  James-Stein  estimates  (see  Chapter  IV),  in 
that  the  estimation  of  rates  for  small  cell  has  improved. 
Also,  an  estimator  ,  designated  the  transformed  scale  cell 
average  (TSCA)  which  is  a  version  corresponding  to  zero 
James-Stein  shrinkage,  has  been  shown  to  be  an  efficient 
estimation  technique,  often  outperforming  all  other  schemes 
examined  here.  The  various  methods  are  quite  competitive, 
and  at  this  time  there  is  no  clear  choice.  We  believe  that 
better  aggregation  methods  need  to  be  applied  prior  to 
attempting  to  choose  among  these  methods. 

E.  ORGANIZATION 

Chapter  II  contains  the  details  of  methodology  and  nota¬ 
tion  necessary  to  the  present  work.  A  brief  summary  of 
James-Stein  estimation  is  presented,  with  emphasis  on  its 
implementation  in  the  present  work. 

Chapter  III  explains  the  limited  translation  extension, 
together  with  the  theoretical  curves  that  help  anticipate 
the  effect  of  this  option. 

Chapter  IV  contains  the  numerical  summaries  and  tabula¬ 
tion  of  the  figures  of  merit  (FOM)  for  the  various 
techniques . 
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Chapter  V  thoroughly  discusses  the  results,  including 
recommendations,  and  lists  additional  areas  needing 
examination. 

The  appendices  document  certain  details  of  interest  to 
the  reader  desiring  a  greater  in-depth  knowledge  of  the 
methods . 


II.  RELATED  ESTIMATION  METHODS 


A.  GENERAL 

This  chapter  describes  the  estimation  method  currently 
in  use  by  the  Marine  Corps,  and  four  other  methods  of 
estimation  pertinent  to  the  present  work. 

B .  BACKGROUND 

As  explained  by  Tucker  [Ref.  1],  the  performance  of  the 
estimation  schemes  are  compared  on  two  scales,  transformed 
and  original.  The  transformed  scale  is  the  range  space  of 
the  Freeman-Tukey  transform  (see  Appendix  B)  that  helps 
stabilize  the  variance  of  the  ordinary  empirical  rates 
assuming  they  are  described  by  the  binomial  model.  On  this 
scale  the  transformed  quantities  are  treated  as  normally 
distributed  random  variables  with  common  variance.  It  is  in 
this  setting  that  the  James-Stein  estimator  is  derived  and 
can  be  expected  to  perform  well.  Because  the  rates  are  low 
and  because  many  cells  are  small,  we  cannot  assert  with 
confidence  that  the  rates  on  the  transformed  scale  are 
approximately  normal  with  common  variance.  The  ultimate 
I  value  must  be  judged  in  terms  of  cross-validation,  i.e., 

comparing  the  estimates  with  like  transformed  values  of 
future  actuals.  Following  Tucker  [Ref.  1]  we  have  chosen  the 
first  four  years  of  data  of  the  seven  available  to  estimate 
I  rates,  and  the  last  three  years  for  validation. 

Although  comparisons  on  the  transformed  scale  are 
valuable  for  purposes  of  understanding  the  behavior  of 
shrinkage  estimators,  they  do  not  supplant  the  need  to  study 
I  behavior  on  the  original  scale.  Hence,  the  transformed 

estimators  must  be  inverted  to  estimated  rates  on  the  orig¬ 
inal  scale,  and  then  validated  against  original  scale 
actuals.  The  traditional  chi-square  goodness  -  of  -  fit 
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statistic  was  chosen  to  do  this.  It  is  a  weighted  sum  of 
squares  deviations  measure. 

The  requirement  to  work  on  a  transformed  scale  intro¬ 
duces  some  additional  complications.  To  illustrate  the 
point,  consider  the  following  choice.  Should  the  empirical 
rate  for  each  cell  of  the  four  years  be  transformed  to  the 
new  scale,  or  should  we  first  sum  the  leavers  and  the  inven¬ 
tory  over  the  four  years  in  order  to  produce  a  more  stable 
cell  rate  prior  to  transforming?  It  turns  out  that  if  the 
latter  is  chosen  then  we  have  no  reasonable  way  to  estimate 
the  within-group  variance  on  the  transformed  scale.  Hence 
the  former  is  chosen.  This  done,  the  transformed  quantities 
are  averaged  over  time  to  produce  a  single  figure  for  each 
cell  prior  to  shrinkage.  (This  is  the  TSCA  mentioned  in 
Chapter  I,  Section  B.) 

On  the  other  hand,  the  latter  figure  is  still  useful 
since  it  is  the  maximum  likelihood  estimator  of  the  cell 
rate  on  the  original  scale.  But  there  is  still  the  question 
of  how  to  use  it,  for  comparison  purposes,  on  the  trans¬ 
formed  scale.  We  have  chosen  to  use  the  four  year  average 
cell  inventory  in  conjunction  with  this  MLE  rate  and  then 
apply  the  arcsine  transformation.  See  equation  B.l. 

To  avoid  confusion,  an  estimate  will  always  be  referred 
to  by  the  name  given  when  initially  calculated.  For  example, 
the  maximum  likelihood  estimate  (MLE)  is  calculated 
initially  on  the  original  scale,  and  still  will  be  called 
the  MLE  when  on  the  transformed  scale.  Also,  the  term 
"maximum  likelihood"  may  be  used  to  refer  to  maximum  likeli¬ 
hood  estimation  in  the  setting  at  hand.  Thus  the  TSCA  ia  a 
set  of  specific  maximum  likelihood  estimators  as  it  refers 
to  the  "IID  normal  with  common  variance"  setting  imposed  on 
the  transformed  scale. 


C.  AGGREGATION 


The  current  Marine  Corps  attrition  rate  estimation  meth¬ 
odology  described  in  Chapter  I  is  an  aggregation  scheme,  and 
as  such  has  several  weaknesses.  First,  a  single  rate  is 
applied  to  all  cells  in  the  aggregate.  This  does  not  take 
into  account  actual  differences  in  cell  inventory  or  losses. 
The  pattern  of  losses  can  differ  greatly  among  cells  in  an 
aggregate  whose  composition  is  arbitrarily  determined. 

Also,  an  aggregation  scheme  has  difficulty  handling 
cells  that  have  had  zero  inventory  for  the  estimation 
period,  i.e.,  cells  whose  MLE  would  be  zero.  The  application 
of  the  aggregate  rate  in  this  case  is  clearly  an 
overestimate . 

Two  constant  rates  are  used  in  the  present  work  for 
comparison  purposes.  The  first  is  the  aggregate  rate  calcu¬ 
lated  on  the  original  scale,  called  hereafter  the  original 
scale  aggregate.  This  rate  is  the  ,  otal  losses  divided  by 
the  total  inventory,  to  be  applied  ..o  all  cells  in  the 
aggregate,  and  is  a  single  number.  T'  r  comparison  on  the 
transformed  scale,  this  rate  is  mapped  into  the  transformed 
space  using  the  arcsine  transformation 

=  (Ni+0.5)^sin'1(2p-l) ,  i  =  l, - K  (2.1) 

where  p  is  the  aforementioned  single  rate,  and  N1  are  the 
average  cell  inventories  over  time.  This  results  in 
different  transformed  cell  means  (transformed  rates)  because 
of  differing  cell  inventories.  Generally  the  subscript  i 
indexes  a  combination  of  LOS  and  OF. 

The  second  rate  is  calculated  on  the  transformed  scale, 
and  is  called  hereafter  the  transform  aggregate.  It  is 
computed  by  averaging  the  cell  figures  that  result  from 
applying  the  the  Freeman- Tukey  transformation.  On  the  trans¬ 
formed  scale  this  is  a  single  number  and  in  fact  is 


grand  mean  of  the  TSCA  figures.  When  inverted  onto  the 
original  scale  for  comparison  purposes  using  equation  A.  12 
of  Ai  pendix  A,  the  average  cell  inventory  is  used,  resulting 
in  d:'  'ferent  rates  for  each  cell. 

D.  MAXIMUM  LIKELIHOOD  ESTIMATION  (MLE ) 

The  present  work  calculates  an  MLE  for  comparison 
purposes.  As  stated  in  Section  B  above,  this  rate  is  the 
total  leavers  (over  time)  divided  by  the  total  inventory 
(over  time),  using  the  four  year  estimation  period.  This 
empirical  rate  is  called  MLE  because  it  would  be  the  maximum 
likelihood  estimate  in  the  setting  of  independent  Bernoulli 
trials.  We  retain  this  terminology  on  the  transformed  scale. 

It  is  well  known  that  this  MLE  is  the  best  unbiased 
estimator  if  the  Bernoulli  setting  is  tenable.  The  problems 
using  the  MLE  here  are  threefold.  First,  the  smaller  the 
number  of  cell  trials  the  greater  the  variability  in  the 
estimation.  Thus  while  the  estimate  is  unbiased,  the  range 
of  values  the  estimate  can  easily  assume  is  large.  A  stable 
estimate  cannot  be  made. 

Second,  this  MLE  assumes  each  data  set  is  drawn  from 
identical  populations.  Service  retention  is  greatly  affected 
by  changing  economic,  political,  and  social  forces.  These 
forces  are  not  constant,  and  over  a  period  of  three  to  six 
years  the  behavior  of  a  cell  can  change  radically.  This 
introduces  the  yearly  update  problem.  The  requirement  to 
drop  old  data  as  new  data  becomes  available  keeps  the  cell 
trials  low,  and  the  variability  of  the  estimate  high.  In 
short,  we  are  not  yet  equipped  to  build  an  estimation  scheme 
based  upon  manpower  flow  model  structural  conditions. 

Third,  manpower  planning  cannot  focus  on  an  individual 
OF/LOS/grade  category.  Since  Marine  Corps  officers  approxi¬ 
mate  a  hierarchical  system,  the  requirement  to  promote  an 
officer  to  fill  a  projected  loss  of,  for  example,  an 
infantry  Lieutenant  Colonel  will  "ripple”  down  to  a 
requirement  to  recruit  an  infantry  Second  Lieutenant. 


E.  TRANSFORMED  SCALE  CELL  AVERAGE  (TSCA) 

Tucker  [Ref.  1]  and  the  present  work  calculate  a  TSCA 
rate.  This  rate  can  be  viewed  as  a  maximum  likelihood  esti¬ 
mator  calculated  on  the  transformed  scale,  using  transformed 
inventory  and  loss  data.  If  the  "normal  with  common  vari¬ 
ance"  model  were  firmly  defensible  on  the  transformed  scale, 
the  TSCA  would  provide  the  best  linear  unbiased  estimators 
of  the  individual  cell  means.  The  method  is  accorded  sepa¬ 
rate  treatment  because  of  the  excellent  results,  especially 
in  the  near  term  (one  year)  validations. 

Stein  [Ref.  3]  in  1955  examined  the  performance  of  this 
maximum  likelihood  estimator  in  predicting  cell  values.  He 
established  that  if  the  number  of  cells  is  at  least  three, 
then  maximum  likelihood  estimation  can  be  improved  in  an 
overall  sense.  The  criterion  he  used  was  the  global  loss, 


L(0  ,  a)  =  p(9i’ai) 


i  =  1,  .  .  .  ,k 


(2.2) 


where  0  is  the  array  of  unobservable  true  cell  values  and  a 
is  the  array  of  predicted  cell  values.  This  global  loss  is 
the  sum  across  all  cells  if  the  individual  loss  is 


=  (01-ai)' 


(2.3) 


where  0^  and  a^  are  the  appropriate  values  for  the  i1"*1  cell 
in  the  aggregate. 

F.  JAMES- STEIN  ESTIMATION 

James  and  Stein  [Ref.  4]  developed  an  estimator,  called 
the  James-Stein  estimator,  which  reduces  the  expected  value 
of  the  global  loss,  when  compared  to  the  cell  means.  The 
expected  value  is  called  the  risk  R: 


=  E  [L(0 ,  a )  ]  =  pt^-a^  ] 


(2.4) 


Two  assumptions  were  made: 

(1)  the  cell  values  are  normally  distributed,  and 

(2)  the  within  cells  variance  is  constant. 

The  basic  idea  of  James-Stein  estimation  is  the  farther 
the  cell  mean  is  from  the  overall  mean,  the  greater  is  the 
size  of  the  residual  error.  Note  that  the  cell  mean  is  the 
TSCA,  and  the  overall  mean  is  the  transform  aggregate  ,  or 
grand  mean.  All  means  are  moved,  or  shrunk,  toward  the 
grand  mean.  The  amount  of  shrinkage  is  proportional  to  the 
absolute  distance  from  the  grand  mean,  i.e.,  the  greater  the 
absolute  distance,  the  greater  is  the  shrinkage. 

There  are,  however,  problems  with  this  method.  First,  a 
natural  objection  is  that  some  cell  attrition  rates  may  be 
far  from  the  grand  mean  simply  because  the  long  term  attri¬ 
tion  from  these  cells  differ  greatly  from  the  majority  of 
cells  in  the  same  aggregate.  To  shrink  the  attrition  rates 
of  these  cells  toward  the  grand  mean  may  be  erroneous. 

In  dealing  with  a  sample  in  the  original  scale,  problems 
occur  when  there  is  no  cell  loss  for  the  entire  estimation 
period,  i.e.,  when  the  MLE  is  zero.  Tucker  [Ref.  1]  handled 
this  by  omitting  such  cells  from  all  comparisons.  In  the 
present  work  these  zero  loss  cells  are  retained  in  an  effort 
to  view  the  effect  of  the  various  schemes  on  the  small  cell 
rate  estimations. 

Third,  Appendix  A  demonstrates  the  Freeman-Tukey  trans¬ 
formation  does  not  normalize  the  cell  means  or  stabilize  the 
variance  when  the  inventory  or  loss  rates  are  low.  Since 
normality  with  common  variance  is  the  basic  assumption  of 
the  James-Stein  scheme,  the  reliability  of  the  results  must 
be  questioned. 

See  Appendix  A  for  the  James-Stein  estimation  algorithm 
as  used  by  Tucker  [Ref.  1]. 

See  James  and  Stein  [Ref.  4]  and  Tucker  [Ref.  1]  for 
details  on  James-Stein  estimation. 


G.  ROBUST  PARAMETRIC  EMPIRICAL  BAYES  ESTIMATION 

An  alternative  method  of  analysis,  unrelated  to  the 
schemes  investigated  here,  is  the  robust  parametric  empir¬ 
ical  Bayes  (RPEB)  model,  suggested  by  D.  P.  Gaver.  While 
relatively  new,  the  method  has  shown  promise  in  the  settings 
to  which  it  has  been  applied. 

Because  of  time  constraints,  this  model  was  not  imple¬ 
mented.  The  procedure  is  a  significant  departure  from  the 
present  work,  but  it  may  offer  important  benefits  to  small 
cell  estimation.  Appendix  F  is  a  brief  description  of  the 
technique . 


III.  LIMITED  TRANSLATION  JAMES- STEIN  ESTIMATION 


A.  GENERAL 

This  chapter  discusses  the  limited  translation  model, 
and  the  validation  approaches  taken  in  this  paper. 

As  stated  in  Chapter  I,  it  is  intuitively  unsettling  to 
shrink  all  empirical  cell  rates  toward  the  grand  mean  by  the 
same  affine  translation.  Also,  one  must  question  whether  the 
risk  could  be  further  reduced  from  that  of  the  James-Stein 
estimator.  Two  articles  published  by  Efron  and  Morris 
[Refs.  5,6]  present  such  a  method:  limited  translation  of 
the  James-Stein  estimator. 

To  compromise  between  James-Stein  and  TSCA  estimation, 
and  to  limit  the  translation  of  extreme  values,  an  interval 
[-C,C]  centered  about  the  grand  mean  is  established.  Inside 
this  interval  all  rates  are  translated  using  full 
James-Stein  shrinkage.  Outside  this  interval  the  amount  of 
shrinkage  is  reduced  the  farther  cell  values  get  from  the 
interval.  The  shrinkage  approaches  zero  in  the  limit. 

To  get  an  intuitive  feel  for  the  differences  between 
James-Stein  and  limited  translation  James-Stein  estimation, 
compare  Figures  2.1  and  3.1.  Figure  2.1  shows  bow  the 
James-Stein  technique  shrinks  all  values  toward  the  grand 
mean.  Figure  3.1  shows  how  limited  translation  estimation 
reduces  the  shrinkage  outside  a  certain  range  of  values 
centered  about  the  grand  mean. 

Theoretically,  limited  translation  estimation,  by 
shrinking  some  cells,  will  slightly  increase  the  global  risk 
over  that  of  the  James-Stein  estimator.  This  increase  is 
acceptable  since  the  individual  cell  risk  of  extreme  inven¬ 
tory  cells  is  decreased.  This  means  the  estimators  of  the 
small  cell  attrition  rates  improve,  usually  significantly, 
at  a  small  cost  to  the  middle  inventory  cells. 
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THE  SHRINKAGE  FACTOR 

1.  The  Limited  Translation  Algorithm 


Conversion  of  the  James-Stein  estimation  algorithm 
to  limit  the  translation  is  straightforward.  From  Appendix 
A,  after  the  attrition  rates  have  been  transformed  using  the 
Freeman-Tukey  transformation,  the  James-Stein  estimator  Pj 
is  calculated, 


pj  =  X  ♦  CjfXi.-X), 


(3.1) 


where 


Cj  =  1  -  (K-3)SSE/ [(K(T-1)+2)SSB] 


(3.2) 


is  the  shrinkage  factor. 

To  modify  Cj  for  limited  translation,  let 


p(u)  -  minimum  (l,d/u'-') 


(3.3) 


where 


u  =  (Xi#-X)  / ( A  *  1 ) 


(3.4) 


Note  that 


A  =  (K(T-1)+2)SSB/(K-3)SSE  -  1 


(3.5) 


and  is  the  variance  of  the  prior  distribution  of  0-  The 
value  of  d  is  chosen  from  the  interval  [0,oc]. 

The  new  shrinkage  factor  is 


CLJ  "  1  '  p(u)[(K-3) SSE/ ((K(T-1)+2)SSB)] 


(3.6) 


and  the  new  estimator  is 


PLJ  =  X  +  CLj(Xi-’X)-  <3-7) 

If  d/u^  is  equal  to  or  greater  than  1,  then  full  shrinkage 
occurs.  However,  if  d/u^  is  less  than  1,  then  full 
shrinkage  does  not  occur. 

2 .  The  Shrinkage  Interval 

The  choice  of  a  d  value  is  important.  The  larger  d 
is,  the  larger  the  interval  [-C,C]  becomes.  If  d  becomes 
infinite,  then  the  estimator  P^j  is  simply  the  James-Stein 
estimator  Pj.  This  would  result  in  no  change  in  the  global 
or  individual  risks.  However,  if  d/u^  is  less  than  one  for 
some  cell  values,  then  those  cells  will  not  fully  shrink, 
and  there  will  be  an  improvement  in  the  small  cell  indi¬ 
vidual  risk.  If  d  equals  zero,  then  the  interval  [-C,C] 
shrinks  to  a  point  at  the  value  of  the  grand  mean,  and  all 
cell  values  shrink  to  the  grand  mean.  In  this  case  the  vari¬ 
ability  of  the  cell  rates  is  zero  Thus,  we  want  some  inter¬ 
mediate  value  of  d.  Appendix  E  discusses  the  theoretical 
implications  of  d,  and  methods  of  choosing  values. 

The  effects  of  d  can  be  directly  observed  in  the 
pattern  of  the  limited  translation  shrinkage  factors.  Within 
[-C,C]  ,  the  cell  shrinkage  factors  will  equal  the 
James-Stein  factor.  The  farther  one  gets  from  the  full 
shrinkage  interval,  the  smaller  the  shrinkage  factors  are. 
Note  the  shrinkage  factor  approaches  zero  in  the  limit. 

Caution  must  be  exercised  when  viewing  the  shrinkage 
factor  pattern  in  this  light.  The  interval  [-C,C]  is  an 
interval  of  cell  values,  in  this  case  cell  means.  Therefore, 
if  the  cell  value  pattern  is  not  apparent,  the  pattern  of 
shrinkage  factors  will  also  not  be  apparent. 


C .  VALIDATION 

The  development  of  the  James- Stein  scheme  took  place  in 
a  restrictive  setting.  It  is  useful  to  state  that  setting  in 
a  form  most  closely  associated  with  our  problem: 

Xij  =  Mi  +  eij>  i=l , . . . ,k  and  j=l,...,n  (3.8) 

with  the  {e  — }  as  IID  N(0,a2)  random  variables.  Since  this 
is  the  same  setting  as  that  of  the  one  way  analysis  of  vari¬ 
ance,  it  is  convenient  to  use  the  quantities  and  notation  of 
ANOVA.  Specifically,  let 


SSL  =  JJtfij-Xi.)*,  V  i,j,  (3.9) 

*  J 

SSB  =  nJ(Xi.-X..)2,  V  i,  (3.10) 

i 

Xi#  =  (l/n)^Xij,  V  i ,  j ,  and  (3.11) 

J 

X..  =  (l/k)JXi#5  Vi.  (3.12) 

i 


Since  we  will  be  using  equations  (3.2),  (7.3),  (7.4), 

and  (7.7)  of  reference  5,  we  record  the  identifications 

'Ni  _ 

a2  - 


SSE/ (2  +  k(n-l))  , 


(3.13) 


(3.15) 


Xc  =  Xc.  for  c= 1 , . . . ,k, 


X  =  X.. ,  and  (3.16) 

p(u)  =  minimum( 1 , d/u^) .  (3.17) 

On  the  other  hand  our  cell  attrition  data  is  fairly 
modeled  with  the  binomial  distribution  [Ref.  1]  such  that 

(1)  y^j  =  number  of  leavers  in  cell  i  during  period  j, 

(2)  N—  =  central  inventory, 

and  the  empirical  cell  attrition  rates  are 

Pi  =  (^ij)/(ZNij)«  Vij.  (3.18) 

J  j 

Because  of  the  large  number  of  small  cells,  the  above  empir¬ 
ical  probabilities  are  unstable  and  it  should  be  possible  to 
improve  the  stability  by  shrinking  them. 

If  we  were  to  treat  the  P-^'s  as  the  in  the  ANOVA 

setting,  we  encounter  some  flaws.  Namely,  the  common  vari¬ 
ance  and  normality  assumptions  are  severly  compromised,  and 
there  is  no  obvious  way  to  estimate  <j2 .  The  same  flaws  are 
present  if  we  back  off  and  use  the  cell  empirical  rates, 

yij/Njj  •  (3.19) 

in  the  role  of  X^j . 

In  hopes  of  giving  relief  to  this  problem  we  (see  refer¬ 
ence  1)  have  chosen  to  use  the  Freeman-Tukey  variance 
stabilization  transform 

Zjj  =  0.5(Nij+0.5)%(sin-1[2(yij)/(nij+l)-l] 

+  sin-1[2(yij-l)/(nij+l)-l]) 
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(3.20) 


in  the  role  of  the  j .  For  the  small  cells  we  still  have 
the  flaws  but  not  to  as  great  an  extent  as  before.  See 
Appendix  B. 

This  leads  us  to  the  perform  some  validation  computa¬ 
tions  using  the  transformed  scale,  i.e.,  the  sum  of  squared 
error  risk  function  will  be  computed  with  the  error  input 
given  by  the  difference  between  the  z^t  (for  tea  cross 
validation  year)  and  the  James-Stein  shrunk  value  z^ 
computed  from  (7.7)  of  reference  5,  i.e., 

z±  =  z..  +  C(zi.-z..)  (3.21) 

and 

C  =  (K-3)SSE/(2+K(n-l))SSB.  (3.22) 

These  risks  will  serve  to  tell  us  whether  the  shrinkage 
technique  is  behaving  as  expected  in  spite  of  the  rough 
treatment  given  the  assumptions  in  the  setting. 

It  is  also  important  to  make  validation  computations  in 
the  original  scale.  Even  though  the  risks  in  the  transformed 
scale  look  attractive,  the  attrition  rate  estimates  must  be 
converted  back  to  the  original  scale.  For  this  purpose  we 
propose  the  chi-square  statistic 

£(yit-NitPi)2/NijPi(l-pi),  V  i  (3.23) 

i 

where  y^t  and  are  the  leavers  and  central  inventory 

counts  for  the  iLn  cell  in  the  tcn  validation  period,  and 

p^  =  . 5 [1- sin( z^f (N^+ . 5 )^] ,  and  (3.24) 


Ni  =  U/nJpij 

J 


V  i  ,  j 


(  3  25 
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provided  the  argument  of  the  sine  function  belongs  to  [-jr/2, 
k/2].  Outside  this  range,  we  use  =  0  or  1  according  to 
whether  the  argument  is  below  ~n/2  or  above  n/2,  respec¬ 
tively.  The  above  value,  p^ ,  is  what  we  will  call  the 
James-Stein  attrition  rate  generator. 

Above  we  have  described  two  validation  risk  calcula¬ 
tions,  one  in  the  transformed  scale  and  one  in  the  original 
scale.  Let  us  now  address  the  question  of  "To  what  are  these 
risks  to  be  compared?" .  The  general  answer  is  to  make  like 
calculations  for  the  other  estimation  schemes:  original  and 
transformed  scale  aggregates,  TSCA,  maximum  likelihood,  and 
the  limited  translation.  To  be  specific,  we  must  address 
some  details.  The  first  four  will  be  taken  up  here  and  the 
limited  translation  modification  will  be  discussed  later. 

We  will  discuss  the  easiest  cases  first:  aggregate  and 
maximum  likelihood  in  the  original  scale.  Let  i  index  all 
cells  in  the  aggregate.  We  define  an  indicator  variable 


D, 


then 


1 

0 


if 


Ni>0 


Ni  =  0 


(3.26) 


,  V  i 


k  =  Pi 

i 


V  i. 


(3.27) 


k  is  thus  the  number  of  cells  in  the  aggregate  with  non- zero 
inventory . 

Now,  if  we  define  another  indicator  variable 


DV  = 


1  N^>0  and  P^gCi^O  or  1  (3.28) 

if  ,  V  i 

0  N^  =  0  or  Pmle^1)^0  or  ^ 


then 


k’  =  P’i 

i 


V  i 


(3.29) 
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k’  is  thus  the  number  of  cells  in  the  aggregate  with 
non- zero  inventory  and  the  MLE  not  equaling  zero  or  one. 


The  aggregate  attrition  rate  is 

p  =  (SZyij)/(ISNij)’  V  ij.  (3.30) 

i  J  '  j 

and  this  value  is  inserted  in  place  of  in  the  chi-square 
statistic.  For  the  maximum  likelihood  method  we  use  instead 

pMLE ( 1 )  =  (Iyij)/(ZNij)>  v  i’J-  (3.31) 

j  j 

In  both  cases  we  must  remove  all  terms  in  the  sum  for  which 
is  zero.  This  has  the  effect  of  reducing  k.  Also,  since 
we  may  have  some  PMLg(i)  =  0  or  1,  these  terms  are  removed, 
leading  to  k'  risk  terms  in  the  sum.  It  is  recommended  that 
such  results  be  multiplied  by  k/k'  in  order  to  provide  a 
fair  comparison  with  methods  that  provide  positive  non-unity 
estimators  for  all  cells. 

Turning  to  the  risk  calculation  in  the  transformed 
scale,  there  are  some  alternative  ways  to  decide  what  we 
call  an  aggregate  or  a  TSCA  estimator.  The  question  arises 
because  to  each  cell  there  is  associated  two  numbers,  y  — 
and  . 

If  performance  in  the  transformed  scale  were  the  only 
concern,  we  would  simply  use  z..  for  the  transformed  aggre¬ 
gate  and  z^,  for  the  TSCA  estimator.  This  in  fact  was  done 
by  Tucker  [Ref.  1] .  But  because  of  the  varying  inventory, 
these  values,  z..  and  z^, ,  do  not  correspond  to  p  and 
PMLE (•*■)•  For  purposes  of  the  present  study  it  was  decided 
to  cast  these  two  values  into  the  transformed  scale: 

za(i)  =  0.5(Ni>0.5)is{sin'1[2(Nip)/(Ni  +  l)-l]  (3.32) 

+  sin"1[2(Nip+l)/(Ni+l)-l] } 


where  is  the  average  inventory  for  the  iun  cell.  This 
choice  leads  to  the  original  scale  aggregation  in  the  trans¬ 
formed  scale  varying  form  cell  to  celi.  It  is  fair  to  do 
this  because  the  varying  inventory  allowance  is  a  part  of 
the  process  that  helps  stabilize  the  value  in  the 
transformed  scale. 

For  the  same  reason,  and  writing  p^  for  P^LE (  i ) »  we  use 

zmle(i)  -  0.5(N-0.5)%{sin'1[2(Nipi)/(Ni+l)-l]  (3.33) 

♦  sin' 1[2(Nipi+l)/ (Ni+1)-1] } 

for  the  maximum  likelihood  estimator  in  the  transformed 


IV.  RESULTS 


A.  GENERAL 

This  chapter  displays  various  data  results  for  the 
estimation  methods  compared  in  this  paper. 

B.  SHRINKAGE  FACTORS 

The  shrinkage  factors  for  the  James-Stein  (JS)  and  the 
limited  translation  James-Stein  (LTJS)  methods  are 
presented,  by  aggregate,  in  Table  1  thru  Table  3.  An  aggre¬ 
gate  is  a  specific  OF  group  and  grade,  e.g.,  the  aviation 
First  Lieutenants,  or  the  combat  support  Lieutenant 
Colonels.  One  James-Stein  factor  is  calculated  for  an  entire 
aggregate,  while  the  limited  translation  technique  assigns  a 
factor  to  each  cell  based  on  the  inventory  and  the  chosen 
value  of  d.  If  the  distributional  assumptions  are  reason¬ 
able,  a  properly  selected  d  will  force  the  middle  inventory 
cells  to  have  equal  James-Stein  and  limited  translation 
factors,  with  ..he  limited  translation  factors  getting 
smaller  the  farther  one  is  from  this  middle  inventory  range. 
As  discussed  in  Appendix  B,  these  assumptions  were  not  met. 
Use  of  a  graph  of  relative  savings  loss  versus  d,  for 
different  validation  years,  yielded  inconclusive  results, 
except  for  the  aviation  First  Lieutenant  aggregate.  See 
Appendix  E.  Since  a  point  of  emphasis  of  this  study  is  small 
cell  estimation,  d  values  were  chosen  to  force  the  limited 
translation  shrinkage  factors  to  act  in  accordance  with 
theoretical  patterns,  as  discussed  in  Chapter  III  and  refer¬ 
ence  5.  The  purpose  was  to  evaluate  the  resulting  small 
cell  risk  values.  See  Table  4. 

The  shrinkage  factors  exhibit  interesting  behavior. 
Within  the  interval  [-C,C]  the  limited  translation  shrinkage 
equals  James-Stein.  Once  outside  the  interval,  the  limited 


translation  factor  is  reduced.  For  example,  consider  Table 
2.  The  code  7  (Engineers)  Lieutenant  Colonels  have  a 
James-Stein  shrinkage  of  .1720.  All  limited  translation 
factors  are  less  than  or  equal  to  this.  Those  that  are  equal 
have  cell  means  (on  the  transformed  scale)  within  [-C,C]. 
Note  that  an  isolated  cell  can  fall  within  this  interval, 
such  as  the  factors  for  LOS's  of  15,  17,  20,  and  29.  For 
those  cell  with  means  outside  [-C,C],  the  respective  limited 
translation  shrinkage  is  less  than  .1720.  See  Chapter  III, 
Section  B,  Subsection  2.  This  behavior  is  non-monotone  with 
respect  to  LOS,  and  upon  reflection,  should  be  anticipated. 
The  pattern  of  officer  losses  is  tied  to  contract  expira¬ 
tions,  reduced  promotional  expectations,  and  retirement 
options . 

C.  FIGURES  OF  MERIT 

Table  5  thru  Table  7  display  the  figures  of  merit  for 
the  six  estimation  schemes  in  transformed  and  original 
space,  as  expected,  both  aggregate  estimation  methods  have 
uniformly  higher  risks  than  the  other  techniques. 

Several  points  are  worth  noting.  First,  in  the  trans¬ 
formed  scale,  limited  translation  ranks  first  overall  in 
lowest  risk.  This  is  contrary  to  the  theory  as  developed  by 
Efron  and  Morris  [Refs.  5,6],  As  noted  in  Chapter  3, 
limited  translation  should  lower  the  individual  cell  risks 
of  those  cells  outside  the  interval  [-C,C],  but  at  the  cost 
of  an  increase  in  global  risk. 

When  the  FOM  is  examined  by  grade,  TSCA  is  always  best 
for  First  Lieutenants,  and  James-Stein  is  always  best  for 
Lieutenant  Colonels.  This  effect  appears  to  be  unchanging. 

In  the  original  scale,  rankings  change  over  time.  For 
1981,  TSCA  is  best  overall  and  for  each  grade.  This  is  also 
true  for  1982,  but  to  a  lesser  degree.  In  1983,  limited 
translation  is  the  best  overall  and  for  each  grade. 


TABLE  1 

AVIATION  SHRINKAGE  FACTORS 


CODE  38 


LOS 

1st  LT 

LTCOL 

JS 

.0128 

.0322 

LTJS  0 

.0093 

.0197 

1 

.0093 

.0197 

2 

.0128 

.0197 

3 

.0085 

.0197 

4 

.0072 

.0197 

5 

.0128 

.0197 

6 

.0128 

.0197 

7 

.0128 

.0197 

8 

.0128 

.0197 

9 

.0128 

.0197 

0 

.0128 

.0197 

11 

.0128 

.0197 

12 

.  0128 

.0197 

13 

.0104 

.0197 

14 

.0096 

.0197 

15 

.0093 

.0234 

16 

.0093 

.0263 

17 

.0093 

.0322 

18 

.0093 

.0177 

19 

.0093 

.0127 

21 

.0093 

.0256 

22 

.0093 

.0205 

23 

.0093 

.0322 

24 

.0093 

.0322 

25 

.0093 

.0322 

26 

.  0093 

.  0322 

27 

.0093 

.0274 

28 

.0093 

.0254 

29 

.0093 

.  0250 

30 

.  0093 

.0192 

The  excellent  performance  of  the  TSCA  estimate  is  worth 
noting.  TSCA  may  be  thought  of  as  the  James-Stein  estimator 
with  zero  shrinkage.  Recognizing  that  shrinkage  is  an  esti¬ 
mated  parameter,  we  are  surprised  in  those  cases  for  which 
the  estimated  shrinkage  is  large,  and  TSCA  outperforms 
James-Stein  and  limited  translation.  See  Tables  3  and  6. 

D.  SMALL  CELL  FIGURES  OF  MERIT 

The  risk  associated  with  the  small  cells  was  investi¬ 
gated  to  determine  which  technique  best  predicts  small  cell 


TABLE  2 

COMBAT  SUPPORT  SHRINKAGE  FACTORS 


CODE  7 

CODE 

13 

CODE 

20 

LOS 

1LT 

LTCOL 

1LT 

LTCOL 

1LT 

LTCOL 

JS 

.  1947 

.  1720 

.  1947 

.  1720 

.  1947 

.  1720 

LTJS  0 

.0645 

.  1531 

.0645 

.  1531 

.0645 

.  1531 

1 

.0645 

.  1531 

.0645 

.  1531 

.  0645 

.  1531 

2 

.0424 

.  1531 

.0501 

.  1531 

.  1947 

.  1531 

3 

.0426 

.  1531 

.0459 

.  1531 

.  1001 

.  1531 

4 

.0323 

.  1531 

.0421 

.  1531 

.  1947 

.  1531 

5 

.0454 

.  1531 

.  1383 

.  1531 

.  1947 

.  1531 

6 

.  1947 

.  1531 

.0978 

.  1531 

.  1368 

.  1531 

7 

.  1947 

.  1531 

.0956 

.  1531 

.  1083 

.  1531 

8 

.0656 

.  1531 

.  1947 

.  1531 

.0422 

.  1531 

9 

.  1947 

.  1531 

.  1116 

.  1531 

.0528 

.  1531 

10 

.  1947 

.  1531 

.0404 

.  1531 

.  1539 

.  1531 

11 

.  1947 

.  1531 

.0551 

.  1531 

.  1947 

.  1531 

12 

.  1947 

.  1531 

.  1044 

.  1531 

.  1947 

.  1531 

13 

.  1947 

.  1531 

.  1947 

.  1531 

.  1947 

.  1531 

14 

.  1387 

.  1531 

.  1743 

.  1531 

.  1947 

.  1531 

15 

.  1556 

.  1720 

.0783 

.  1531 

.  1190 

.  1531 

16 

.  1947 

.  1531 

.0645 

.  1720 

.  1947 

.  1720 

17 

.  1947 

.  1720 

.0645 

.  1720 

.  1947 

.  1720 

18 

.  1947 

.0898 

.0645 

.0997 

.  1947 

.  1720 

19 

.  1947 

.0718 

.0833 

.0645 

.  1947 

.  1720 

20 

.  1947 

.  1720 

.0833 

.  1353 

.0965 

.  1720 

21 

.  1947 

.0919 

.0645 

.  1383 

.  1947 

.  1720 

22 

.  1947 

.  1447 

.0645 

.  1600 

.  1277 

.  1720 

23 

.  1947 

.  1575 

.0645 

.  1671 

.  1539 

.  1720 

24 

.  1084 

.  1720 

.0645 

.  1720 

.0980 

.  1531 

25 

.0833 

.  1720 

.0645 

.  1720 

.0995 

.  1531 

26 

.0833 

.  1720 

.0645 

.  1720 

.0995 

.  1531 

27 

.0645 

.  1061 

.0645 

.  1720 

.  0974 

.  1531 

28 

.0645 

.  1418 

.0645 

.  1720 

.0727 

.  1531 

29 

.0645 

.  1720 

.0645 

.  1720 

.0727 

.  1720 

30 

.0645 

.  1555 

.  0645 

.  1720 

.  1077 

.  1253 

attrition  since  the  global  figures  of  merit  in  Tables  4  thru 
6  may  be  disguising  what  is  happening  in  the  small  cells. 
The  small  cell  figures  of  merit  tables  record  the  small  cell 
contribution  to  the  global  risk.  Table  8  thru  Table  13 
display  these  results.  Two  average  inventory  ranges  were 
examined:  zero  to  five,  and  six  to  ten. 

For  the  range  zero  to  five,  in  the  transformed  scale, 
the  results  are  varied.  The  TSCA  method  ranks  first  for 
1981,  limited  translation  first  for  1982,  and  James-Stein 


TABLE  3 

GROUND  COMBAT  SHRINKAGE  FACTORS 


CODE  3 

CODE  5 

CODE 

10 

LOS 

1LT 

LTCOL 

1LT 

LTCOL 

1LT 

LTCOL 

JS 

.0796 

.0392 

.0796 

.0392 

.0796 

.0392 

LTJS  0 

.0480 

.0241 

.  0480 

.0241 

.0480 

.0241 

1 

.0507 

.0241 

.0480 

.0241 

.0480 

.0241 

2 

.0343 

.0241 

.0441 

.0241 

.0796 

.0241 

3 

.0268 

.0241 

.0527 

.0241 

.0796 

.0241 

4 

.0313 

.0241 

.0744 

.  0241 

.0796 

.0241 

5 

.0281 

.0241 

.0641 

.0241 

.0796 

.0241 

6 

.  0796 

.0241 

.0796 

.0241 

.0796 

.0241 

7 

.  0796 

.0241 

.0796 

.0241 

.0732 

.0241 

8 

.0796 

.0241 

.0796 

.0241 

.0796 

.  0241 

9 

.  0796 

.0241 

.0796 

.0241 

.  0796 

.0241 

10 

.0796 

.0241 

.0796 

.0241 

.0796 

.0241 

11 

.0796 

.0241 

.0796 

.0241 

.  0796 

.0241 

12 

.0796 

.0241 

.0796 

.0241 

.  0753 

.0241 

13 

.0796 

.0241 

.0796 

.0241 

.0631 

.0241 

14 

.0796 

.0241 

.0657 

.0241 

.  0507 

.0241 

15 

.  0663 

.0241 

.0537 

.0241 

.0480 

.0241 

16 

.  0654 

.0392 

.  0537 

.  0368 

.  0480 

.0275 

17 

.  0480 

.0217 

.  0480 

.0392 

.  0480 

.0392 

18 

.  0480 

.0092 

.0480 

.0191 

.0480 

.0392 

19 

.0480 

.0082 

.0480 

.0208 

.0480 

.0392 

20 

.0480 

.0147 

.0480 

.0392 

.0480 

.0392  i 

21 

.0480 

.0118 

.0480 

.0291 

.0480 

.0392 

22 

.0480 

.0165 

.0480 

.0392 

.0480 

.0392 

23 

.0480 

.0215 

.0480 

.0392 

.0480 

.0392 

24 

.0480 

.0392 

.0480 

.0386 

.0480 

.0383 

25 

.0480 

.0392 

.0480 

.0392 

.0480 

.0271 

26 

.0480 

.0392 

.0480 

.0233 

.0480 

.0389 

27 

.0480 

.0392 

.0480 

.0221 

.0480 

.0309 

28 

.  0480 

.  0392 

.  0480 

.0267 

.  0480 

.0323 

29 

.0480 

.0392 

.0480 

.  0213 

.  0480 

.0286 

30 

.  0480 

.  0273 

.  0480 

.  0204 

.0480 

.0283  ' 

first  for  1983.  In  the  original  scale  TSCA  is  first  for 
1981,  James-Stein  first  for  1982,  and  limited  translation 
first  for  1983.  The  most  surprising  result  is  that  no  clear 
pattern  has  emerged,  other  than  the  consistently  poor 
performance  of  the  two  aggregate  estimation  schemes.  The 
limited  translation  technique  has  not  resulted  in  a  uniform 
lowering  of  the  small  cell  risk.  In  fact,  in  about  one-half 
of  the  cases  James-Stein  estimation  resulted  in  a  lower 
small  cell  risk  than  did  limited  translation. 


TABLE  4 

D  VALUES  USED  IN  STUDY  AGGREGATES 


1st  LT 

LTCOL 

Aviation 

.8 

.4 

Combat  Support 

.2 

.4 

Ground  Combat 

.5 

.  3 

Table  14  thru  Table  19  list  the  small  cell  percentages 
of  total  risk.  These  results  are  useful  in  determining  the 
actual  contribution  of  the  small  cells  tc  the  total  risk  for 
the  specific  aggregate,  estimation  scheme,  and  year. 

E.  ATTRITION  RATES 

Since  the  ultimate  purpose  of  this  study  is  to  produce 
attrition  rates,  the  rates  generated  by  the  six  estimation 
methods  are  presented  in  Table  20  thru  Table  33.  The 
differing  results  of  each  method  are  apparent.  Of  interest 
is  the  relative  agreement  among  MLE,TSCA,  James-Stein,  and 
limited  translation,  when  compared  to  the  two  aggregate 
methods . 

Also,  we  are  comforted  '•’■'at  the  attrition  rate  patterns 
follow  both  the  pattern  of  the  raw  attrition,  and  experi¬ 
ence.  Experience  tells  us  that,  for  First  Lieutenants,  no 
attrition  occurs  within  LOS's  of  zero  or  one  because  there 
are  no  First  Lieutenants  with  such  LOS’s.  A  peak  in  attri¬ 
tion  should  occur  around  the  eight  year  point,  because  most 
officers  have  been  promoted  to  Captain,  and  those  remaining  I 

resign  to  pursue  civilian  careers.  Another  peak  should  occur 
at  twenty  years  since  at  this  point  First  Lieutenants  who 
have  enlisted  time  start  to  retire. 
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Tucker  [Ref.  1]  demonstrated  this  loss  pattern  in  a 
series  of  graphs  (Figures  A.l  thru  A.  10  in  reference  1) 
which  display  raw  loss  rates  for  selected  grades  and  OF '  s . 
These  graphs  show  an  increase  in  attrition  rate  peaking  at 
eight  years  for  First  lieutenants,  and  again  at  twenty 
years . 

The  rates  for  TSCA,  James-Stein,  and  limited  translation 
follow  this  pattern.  MLE  does  also,  in  general.  However,  MLE 
will  predict  zero  attrition  if  the  cells  have  been  empty  for 
the  entire  estimation  period,  as  it  does  in  Table  22.  The 
predicted  rates  for  LOS  13  thru  19  are  zero,  while  the  above 
three  schemes  estimate  rates  ranging  from  .07  to  .09. 

The  aggregate  methods  do  not  display  this  pattern.  In 
fact,  the  translated  scale  aggregate  estimate  is  at  a  local 
minimum  when  the  above  schemes  reach  local  maximums .  The 
original  scale  aggregate  ,  of  course,  remain  constant  except 
for  those  cells  forced  to  zero  by  zero  cell  inventories. 


TABLE  5 

AVIATION  FIGURES  OF  MERIT 

1981  1982  1983 

TRANSFORMED  FOM 


1st  LT 


AGG  ORIG 

6.405 

9  .  645 

11.393 

AGG  TRANS 

197 . 621 

202 . 022 

208 . 808 

MLE 

3.914 

9 . 981 

10.420 

TSCA 

3.461 

9.574 

10.042 

JS 

3 .678 

9 .768 

10.318 

LTJS 

3  .  642 

9 .764 

10.279 

LTCOL 

AGG  ORIG 

9.957 

17 . 997 

15.488 

AGG  TRANS 

25 . 093 

29 . 506 

29.310 

MLE 

4.366 

9 . 058 

8.210 

TSCA 

5  .777 

10.967 

10.394 

JS 

5  .737 

10.911 

10.355 

LTJS 

5  .  744 

10 . 901 

10.340 

ORIGINAL  FOM 

1st  LT 

AGG  ORIG 

31.499 

33 . 146 

57 . 504 

AGG  TRANS 

310.625 

286 . 068 

566 . 903 

MLE 

24.455 

48 . 375 

57.475 

TSCA 

19.333 

34.678 

49. 105 

JS 

22 . 094 

39 . 596 

51.040 

LTJS 

21.081 

40.590 

50.690 

LTCOL 

AGG  ORIG 

74 . 697 

110.793 

50.456 

AGG  TRANS 

19276.410 

248 . 171 

4894 . 164 

MLE 

34.473 

56 . 304 

22.  116 

TSCA 

37.9  74 

51.410 

29.613 

JS 

43 . 187 

54 . 901 

33.736 

LTJS 

43.802 

57  .  964 

35.619 

MAXIMUM  LIKELIHOOD 

COMPARISON  FACTOR 

k 

k' 

k/k’ 

1st  LT  13 

11 

1.1818 

LTCOL  16 

13 

1 . 2308 

TABLE  6 

COMBAT  SUPPORT  FIGURES  OF  MERIT 

1981  1982  1983 

TRANSFORMED  FOM 


1st  LT 


AGG  ORIG 

1.528 

1.842 

1.385 

AGG  TRANS 

3  .  308 

3  .  722 

3.238 

MLE 

2.273 

2.843 

2.409 

TSCA 

1.383 

2 . 008 

1.468 

JS 

1.470 

2 . 066 

1.532 

LTJS 

1.428 

2.045 

1.482 

LTCOL 

AGG  ORIG 

1.316 

1.701 

1.755 

AGG  TRANS 

2.601 

2 . 898 

2.981 

MLE 

1.589 

2 . 157 

2.443 

TSCA 

.910 

1.689 

1.752 

JS 

.809 

1.514 

1.556 

LTJS 

.831 

1.560 

1.611 

ORIGINAL  FOM 

1st  LT 

AGG  ORIG 

79.521 

74.294 

90.992 

AGG  TRANS 

157 . 641 

131.044 

27231.380 

MLE 

123.508 

128.483 

104.668 

TSCA 

73.190 

80.021 

122.850 

JS 

77.802 

81.133 

91.139 

LTJS 

80.324 

83.379 

68.079 

LTCOL 

48.780 

AGG  ORIG 

34.631 

104.748 

AGG  TRANS 

48.775 

44.739 

180.760 

MLE 

57.249 

66 . 996 

48.468 

TSCA 

27 .773 

40.413 

87.428 

JS 

27.804 

36 . 758 

54.514 

LTJS 

35.847 

37.122 

42.044 

MAXIMUM  LIKELIHOOD 

COMPARISON  FACTOR 

k 

k' 

k/k’ 

1st  LT  70 

33 

2 . 1212 

LTCOL  40 

23 

1.7391 

TABLE  7 

GROUND  COMBAT  FIGURES  OF  MERIT 


1981 

1982 

1983 

TRANSFORMED 

FOM 

1st  LT 

AGG  ORIG 

2.774 

3 . 609 

5  .  720 

AGG  TRANS 

18.280 

19.064 

22.044 

MLE 

2.693 

4.258 

6.178 

TSCA 

1.957 

3.447 

5.231 

JS 

2.169 

3.629 

5.505 

LTJS 

2.045 

3.482 

5.334 

LTCOL 

AGG  ORIG 

3  .  692 

3.596 

3 . 783 

AGG  TRANS 

13 . 108 

12.209 

12.546 

MLE 

1.252 

2.029 

2.925 

TSCA 

1.598 

2.414 

3 . 333 

JS 

1.584 

2.332 

3.223 

LTJS 

1.567 

2.367 

3.264 

ORIGINAL 

FOM 

1st  LT 

AGG  ORIG 

79.406 

80.340 

99 .979 

AGG  TRANS 

288.193 

255.215 

321.524 

MLE 

111.881 

145.280 

170.424 

TSCA 

69 . 166 

92.226 

106 .611 

JS 

83.893 

98.672 

117 . 320 

LTJS 

84.775 

101.875 

121.432 

LTCOL 

AGG  ORIG 

100.321 

118.816 

111.066 

AGG  TRANS 

591.887 

342.375 

30663 .610 

MLE 

42.304 

43.437 

56.710 

TSCA 

37.521 

49.664 

57.169 

JS 

42.460 

50.712 

61.459 

LTJS 

45.192 

50.940 

55.562 

MAXIMUM  LIKELIHOOD  COMPARISON  FACTOR 

k  k’  k/k' 

1st  LT  45  35  1.2857 

LTCOL  46  35  1.3143 


TABLE  8 


AVIATION  SMALL  CELL  FOM 
(INV  £  5) 


1981 

1982 

1983 

TRANSFORMED 

FOM 

1st  LT 

AGG  ORIG 

0.069 

0.021 

0.017 

AGG  TRANS 

6.612 

4.682 

5.956 

MLE 

0.176 

0.001 

0.083 

TSCA 

0.021 

0.069 

0.0001 

JS 

0.031 

0.054 

0.002 

LTJS 

0.029 

0.056 

0.001 

LTCOL 

AGG  ORIG 

0.434 

1.015 

0.493 

AGG  TRANS 

2.612 

3 . 628 

2.507 

MLE 

0.277 

0.640 

0.136 

TSCA 

0.259 

0.617 

0.328 

JS 

0.265 

0.637 

0.308 

LTJS 

0.266 

0.630 

0.314 

ORIGINAL  FOM 

1st  LT 

AGG  ORIG 

5.277 

0.146 

0 . 504 

AGG  TRANS 

0 

0 

0 

MLE 

0 

0 

0 

TSCA 

0.882 

0.520 

0.520 

JS 

1.432 

0.386 

0.386 

LTJS 

0.436 

0.136 

0.136 

LTCOL 

AGG  ORIG 

6.034 

12.751 

13.535 

AGG  TRANS 

98 . 746 

138.781 

4814.621 

MLE 

5.846 

8.400 

1.587 

TSCA 

3  .  355 

6.  112 

5 . 264 

JS 

5 . 960 

8 . 279 

4.637 

LTJS 

4 . 870 

7 .719 

2 . 737 

TABLE  9 

AVIATION  SMALL  CELL  FOM 
(6  S  INV  5  10) 


1981  1982  1983 

TRANSFORMED  FOM 


1st  LT 

AGG  ORIG 

0. 139 

0.136 

0.0001 

AGG  TRANS 

4.960 

4.960 

3 . 500 

MLE 

0.085 

0.085 

0 . 004 

TSCA 

0.072 

0.072 

0 . 008 

JS 

0.086 

0.086 

0 . 004 

LTJS 

0.086 

0.086 

0.004 

LTCOL 

AGG  ORIG 

1.966 

0.588 

0.146 

AGG  TRANS 

2.777 

0.711 

0.195 

MLE 

1.096 

0.134 

0 . 00001 

TSCA 

1.535 

0.448 

0.001 

JS 

1.541 

0.425 

0 . 0004 

LTJS 

1.541 

0.432 

0 . 0004 

1st  LT 

AGG  ORIG 

1.564 

ORIGINAL  FOM 

1.564 

0.331 

AGG  TRANS 

0 

0 

0 

MLE 

0.692 

0 .692 

0 . 556 

TSCA 

0.459 

0.459 

0.575 

JS 

0.643 

0.643 

0 . 507 

LTJS 

0.594 

0.594 

0.468 

LTCOL 

AGG  ORIG 

36.632 

9.338 

0.332 

AGG  TRANS 

19108.090 

40 . 132 

0 . 726 

MLE 

9.495 

0.821 

0.462 

TSCA 

15 . 302 

3.704 

0 .713 

JS 

12 . 999 

2.586 

0 . 680 

LTJS 

9.943 

1.992 

0.453 

TABLE  10 


COMBAT  SUPPORT  SMALL  CELL  FOM 
(INV  £  5) 


1981 

1982 

1983 

TRANSFORMED 

FOM 

1st  LT 

AGG  ORIG 

0.509 

0.305 

0.473 

AGG  TRANS 

0.673 

0.355 

0.518 

MLE 

0.945 

0 . 929 

1.229 

TSCA 

0.436 

0 . 369 

0.632 

JS 

0.432 

0.326 

0.564 

LTJS 

0.445 

0.360 

0.606 

LTCOL 

AGG  ORIG 

0.791 

1.073 

1.055 

AGG  TRANS 

0.614 

0.612 

0.897 

MLE 

1.234 

1.736 

2.164 

TSCA 

0.753 

1.391 

1.479 

JS 

0.650 

1.159 

1.241 

LTJS 

0.670 

1.217  1.308 

ORIGINAL  FOM 

1st  LT 

AGG  ORIG 

35.470 

26.300 

51.054 

AGG  TRANS 

54.725 

19.755 

27137.010 

MLE 

27.405 

13.489 

24.454 

TSCA 

30.167 

24.187 

82.679 

JS 

30.208 

19.083 

46.552 

LTJS 

36.180 

26.113 

29.901 

LTCOL 

AGG  ORIG 

28.480 

41.359 

102.329 

AGG  TRANS 

35.152 

25.355 

167.969 

MLE 

50.620 

58.087 

43 . 565 

TSCA 

22.688 

32.333 

84.557 

JS 

22 . 687 

27.958 

49 .336 

LTJS 

31 . 774 

30.339 

36.854 

TABLE 

11 

COMBAT 

SUPPORT  SMALL  CELL  FOM 
(6  £  INV  S  10) 

1981 

1982 

1983 

TRANSFORMED  FOM 

1st  LT 

AGG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

0.252 

0.203 

0.155 

0.223 

0.196 

0.221 

0.113 

0.151 

0.102 

0.086 

0.082 

0.086 

0.136 
0.136 
0.040 
0.061 
0 . 057 
0.060 

LTCOL 

AGG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

0.112 

0.353 

0.166 

0.083 

0.068 

0.072 

0.143 

0.494 

0.162 

0.177 

0.215 

0.204 

0.302 
0.616 
0 . 141 
0.252 
0.290 
0.278 

ORIGINAL  FOM 

1st  LT 

AGG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

9.745 

32 . 916 
24.051 
8.198 
8.697 
7.967 

9.120 

9  .  Ill 
15.102 
6.275 
6.188 
5.314 

8.899 
15.834 
6  459 
2.848 
3.580 
2.469 

LTCOL 

AGG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

5.572 

9.841 

4 . 665 
3.173 
3.397 
2.809 

4.976 

11.651 

6.313 

6.486 

5.871 

4.280 

2.419 
12.791 
4.904 
2.871 
5.177 
5  .  190 

TABLE  12 


GROUND  COMBAT  SMALL  CELL  FOM 
(INV  £  5) 


1981 

1982 

1983 

TRANSFORMED 

FOM 

1st  LT 

AGG  ORIG 

0.246 

0.104 

0 . 209 

AGG  TRANS 

1.810 

1.392 

1.683 

MLE 

0.447 

0.365 

0.528 

TSCA 

0.245 

0.161 

0.285 

JS 

0.252 

0.171 

0 . 302 

LTJS 

0.253 

0.170 

0.300 

LTCOL 

AGG  ORIG 

0.762 

0.629 

0.569 

AGG  TRANS 

1.692 

1.703 

1.660 

MLE 

0.396 

0.450 

0.501 

TSCA 

0  .’380 

0.388 

0.459 

JS 

0 .376 

0.388 

0.446 

LTJS 

0.377 

0.385 

0.451 

ORIGINAL 

FOM 

1st  LT 

AGG  ORIG 

9.413 

7  .  352 

16.442 

AGG  TRANS 

44.715 

17.024 

72.074 

MLE 

8.221 

7.687 

20. 185 

TSCA 

10.589 

8.740 

17.879 

JS 

11.901 

3.733 

13.206 

LTJS 

9.775 

3.916 

11.794 

LTCOL 

AGG  ORIG 

40.102 

41.504 

59 . 156 

AGG  TRANS 

481.735 

200.544 

30543.860 

MLE 

24.230 

15.215 

24.940 

TSCA 

13 . 307 

16 . 803 

24.049 

JS 

18.610 

16 .513 

25 . 658 

LTJS 

22 . 270 

16.885 

17 . 772 

TABLE  13 


GROUND  COMBAT  SMALL  CELL  FOM 
(6  S  INV  5  10) 


1981  1982  1983 

TRANSFORMED  FOM 


1st  LT 


AGG  ORIG 

0.097 

0.150 

0.160 

AGG  TRANS 

0.054 

0.115 

0.131 

MLE 

0.011 

0.180 

0.251 

TSCA 

0.010 

0.173 

0.249 

JS 

0.012 

0.163 

0.237 

LTJS 

0.012 

0.163 

0.237 

LTCOL 

1.114 

0.394 

AGG  ORIG 

0 .530 

AGG  TRANS 

0.224 

0.696 

0.204 

MLE 

0.295 

0.462 

0.600 

TSCA 

0.535 

0.799 

0.615 

JS 

0.509 

0.777 

0.578 

LTJS 

0.514 

0.780 

0.588 

ORIGINAL 

FOM 

1st  LT 

AGG  ORIG 

5.659 

3.658 

5.081 

AGG  TRANS 

3.490 

2.767 

4.614 

MLE 

2.809 

5.387 

13.237 

TSCA 

2.310 

3.265 

8.005 

JS 

2.501 

3.970 

8 . 155 

LTJS 

2.664 

3.904 

9.178 

LTCOL 

47.989 

AGG  ORIG 

9.408 

10.095 

AGG  TRANS 

5.452 

43.834 

8  .  123 

MLE 

6.041 

13.789 

13.300 

TSCA 

11.155 

17.113 

10 . 940 

JS 

9.481 

18.416 

11.942 

LTJS 

8.467 

17.962 

13 . 923 

TABLE  14 


AVIATION  SMALL  CELL  PERCENTAGE  OF  TOTAL  FOM 

(INV  :£  5) 


1981 

1982 

TRANSFORMED 

1983 

1st  LT 

AGG  ORIG 

0.011 

0.002 

0.002 

AGG  TRANS 

6.033 

4.023 

5 . 029 

MLE 

3.045 

0.00001 

0.008 

TSCA 

0.006 

0.007 

0 . 00002 

JS 

0.008 

0.005 

0.0002 

LTJS 

0.008 

0.006 

0.0002 

LTCOL 

AGG  ORIG 

0.044 

0.056 

0.032 

AGG  TRANS 

0. 104 

0. 123 

0.086 

MLE 

0.063 

0.071 

0.017 

TSCA 

0.045 

0 . 056 

0.032 

JS 

0.046 

0 . 058 

0.030 

LTJS 

0.046 

0.057 

ORIGINAL 

0 .030 

1st  LT 

AGG  ORIG 

0 . 168 

0.005 

0 . 003 

AGG  TRANS 

0 

0 

0 

MLE 

0 

0 

0 

TSCA 

0.046 

0.014 

0.010 

JS 

0.065 

0.010 

0.008 

LTJS 

0.021 

0.003 

0.003 

LTCOL 

AGG  ORIG 

0.081 

0.115 

0.068 

AGG  TRANS 

0.005 

0.560 

0 . 984 

MLE 

0 . 170 

0. 149 

0 .072 

TSCA 

0.088 

0. 119 

0. 178 

JS 

0. 138 

0.151 

0 . 137 

LTJS 

0.111 

0.116 

0.077 

TABLE  15 


AVIATION  SMALL  CELL  PERCENTAGE  OF  TOTAL  FOM 
(6  £  INV  2  10) 


1981 

1982 

1983 

1st  LT 

TRANSFORMED 

AGG  ORIG 

0.021 

0.014 

0.00001 

AGG  TRANS 

0.025 

0.025 

0.017 

MLE 

0.022 

0.009 

0.0004 

TSCA 

0.021 

0.007 

0.0008 

JS 

0.023 

0.009 

0 . 0004 

LTJS 

0.024 

0.009 

0.0004 

LTCOL 

AGG  ORIG 

0 . 197 

0.033 

0.009 

AGG  TRANS 

0.  Ill 

0.024 

0 . 007 

MLE 

0.251 

0.015 

2E-  6 

TSCA 

0 . 266 

0.051 

0.0001 

JS 

0.269 

0.039 

0.00003 

LTJS 

0.268 

0.040 

0.00003 

1st  LT 

ORIGINAL 

AGG  ORIG 

0.050 

0.047 

0.006 

AGG  TRANS 

0 

0 

0 

MLE 

0.028 

0.014 

0.010 

TSCA 

0.024 

0.013 

0.012 

JS 

0.029 

0.016 

0.010 

LTJS 

0.028 

0.014 

0.009 

LTCOL 

AGG  ORIG 

0.490 

0.084 

0.007 

AGG  TRANS 

0.991 

0.162 

0.0001 

MLE 

0.275 

0.015 

0.021 

TSCA 

0.403 

0.072 

0.024 

JS 

0.301 

0.047 

0.020 

LTJS 

0.227 

0.034 

0.013 

TABLE  16 

COMBAT  SUPPORT  SMALL  CELL  PERCENTAGE  OF  TOTAL  FOM 

(INV  5  5) 


1981 

1982 

TRANSFORMED 

1983 

i  1st  LT  ... 

AGG  ORIG 

0.334 

0. 165 

0.342 

AGG  TRANS 

0.203 

0.095 

0 . 160 

MLE 

0.416 

0 . 327 

1.510 

TSCA 

0.316 

0 . 184 

0.430 

JS 

0.294 

0.158 

0.368 

LTJS 

0.312 

0 . 176 

0.409 

1  LTCOL  .... 

AGG  ORIG 

0.602 

0.631 

0.601 

AGG  TRANS 

0.236 

0 .211 

0 . 301 

MLE 

0 .776 

0.805 

0 . 886 

TSCA 

0.828 

0.823 

0.844 

JS 

0 . 803 

0.726 

0.798 

LTJS 

0.806 

0.782 

ORIGINAL 

0.812 

!  1st  LT 

AGG  ORIG 

0.446 

0.354 

0 .561 

AGG  TRANS 

0.347 

0.151 

0.997 

MLE 

0.222 

0.  105 

0 . 234 

TSCA 

0.412 

0.302 

0 .763 

!  js 

0.388 

0.235 

0.511 

i  LTJS 

0.450 

0.313 

0.439 

1  LTCOL  _ 

AGG  ORIG 

0.822 

0 . 848 

0 . 977 

AGG  TRANS 

0.721 

0.567 

0 . 929 

|  MLE 

0 . 884 

0 .867 

0.899 

TSCA 

0.817 

0 . 800 

0.967 

!  JS 

0.816 

0.761 

0 . 905 

LTJS 

0 .886 

0.817 

0 .877 

TABLE  17 


COMBAT  SUPPORT  SMALL  CELL  PERCENTAGE 

(6  S  INV  £  10) 

1981 


1st  LT 


A GG  ORIG 

0. 165 

AGG  TRANS 

0 . 061 

MLE 

0.068 

TSCA 

0. 161 

JS 

0.134 

LTJS 

0.155 

LTCOL 

AGG  ORIG 

0.085 

AGG  TRANS 

0 . 136 

MLE 

0 . 105 

TSCA 

0.091 

JS 

0.084 

LTJS 

0.087 

1st  LT 

AGG  ORIG 

0. 123 

AGG  TRANS 

0.209 

MLE 

0.195 

TSCA 

0.  112 

JS 

0.112 

LTJS 

0.099 

LTCOL 

AGG  ORIG 

0.161 

AGG  TRANS 

0.202 

MLE 

0.081 

TSCA 

0. 114 

JS 

0. 122 

LTJS 

0 .078 

OF  TOTAL  FOM 


1982 

1983 

TRANSFORMED 

0.061 

0  098 

0.041 

0.042 

0.036 

0.017 

0.043 

0.041 

0 . 040 

0.037 

0.042 

0 . 041 

0 . 084 

0.172 

0.017 

0 . 207 

0.075 

0 . 058 

0.105 

0.144 

0 . 142 

0 . 186 

0. 131 

0.173 

ORIGINAL 

0. 123 

0.098 

0.070 

0 . 0006 

0.118 

0  062 

0.078 

0.023 

0.076 

0.039 

0.064 

0.036 

0.102 

0.023 

0.260 

0.071 

0.094 

0.101 

0.161 

0.033 

0.160 

0 . 095 

0.  115 

0.123 

TABLE  18 


GROUND  COMBAT  SMALL  CELL  PERCENTAGE  OF  TOTAL  FOM 

(INV  5  5) 


1981 

1982 

1983 

TRANSFORMED 

1st  LT 

AGG  ORIG 

0 . 089 

0.029 

0 . 037 

AGG  TRANS 

0.099 

0.073 

0 . 076 

MLE 

0. 166 

0.086 

0.085 

TSCA 

0.125 

0.047 

0 . 054 

JS 

0 . 116 

0.047 

0 . 055 

LTJS 

0 . 124 

0.049 

0 . 050 

LTCOL 

AGG  ORIG 

0.206 

0.175 

0.150 

AGG  TRANS 

0. 129 

0.  139 

0. 132 

MLE 

0.316 

0 . 222 

0.171 

TSCA 

0.248 

0 . 161 

0  .  138 

JS 

0.248 

0 . 166 

0 . 138 

LTJS 

0.240 

0 . 163 

0.  138 

ORIGINAL 

1st  LT 

AGG  ORIG 

0.119 

0.092 

0 . 164 

AGG  TRANS 

0  .  115 

0.067 

0.224 

MLE 

0.073 

0 . 053 

0.118 

TSCA 

0.153 

0.075 

0. 168 

JS 

0.  143 

0.038 

0.  113 

LTJS 

0.117 

0.038 

0.097 

LTCOL 

AGG  ORIG 

0 . 400 

0 . 349 

0.533 

AGG  TRANS 

0.814 

0 . 586 

0 . 996 

MLE 

0.574 

0  .  350 

0.440 

TSCA 

0 . 345 

0.338 

0.421 

JS 

0 .438 

0.326 

0.417 

LTJS 

0.493 

0.331 

0.320 

TABLE  19 


GROUND  COMBAT  SMALL  CELL  PERCENTAGE  OF  TOTAL  FOM 

(6  S  INV  <;  10) 


1981 

1982 

TRANSFORMED 

1983 

1st  LT 

AGG  ORIG 

0 .035 

0.042 

0.028 

AGG  TRANS 

0.003 

0 . 006 

0.006 

MLE 

0.004 

0.042 

0.041 

TSCA 

0 . 005 

0.050 

0.048 

JS 

0.005 

0.045 

0.043 

LTJS 

0 . 006 

0.047 

0.044 

LTCOL 

AGG  ORIG 

0.144 

0.310 

0.104 

AGG  TRANS 

0 .017 

0.057 

0.016 

MLE 

0.236 

0.223 

0.205 

TSCA 

0.335 

0.331 

0. 185 

JS 

0.322 

0.333 

0 . 179 

LTJS 

0 . 328 

0.330 

ORIGINAL 

0.180 

1st  LT 

AGG  ORIG 

0.071 

0.046 

0.051 

AGG  TRANS 

0.012 

0.011 

0 .014 

MLE 

0.025 

0.037 

0.078 

TSCA 

0.003 

0.035 

0.075 

JS 

0.030 

0.040 

0.070 

LTJS 

0.031 

0.038 

0.076 

LTCOL 

AGG  ORIG 

0.094 

0.404 

0.091 

AGG  TRANS 

0.009 

0.128 

0.0003 

MLE 

0.143 

0.317 

0.235 

TSCA 

0 . 297 

0.345 

0.191 

JS 

0 . 223 

0.363 

0.194 

LTJS 

0. 187 

0.353 

0.251 

TABLE 

20 

AVIATION 

ATTRITION 

RATES 

FOR  1ST  LTS 

.OS 

AGG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

2 

.0485 

.  1079 

.0150 

.0160  . 

0167 

.0167 

3 

.0485 

.2218 

.0265 

.0277 

0291 

.0286 

4 

.0485 

.2466 

.0304 

.0310 

0325 

.0318 

5 

.0485 

.2155 

.0954 

.0901 

0914 

.0914 

6 

.0485 

.  1367 

.0488 

.0490 

0498 

.0498 

7 

.0485 

.0687 

.  0600 

.0609 

0610 

.0610 

8 

.0485 

.0037 

.0724 

.0654 

0641 

.0641 

9 

.0485 

.0039 

.0443 

.0474 

0465 

.0465 

10 

.0485 

.0043 

.0659 

.0704 

0691 

.0691 

11 

.0485 

.0255 

.0861 

.0871 

0838 

.0838 

12 

.0485 

0 

.0727 

.0812 

0723 

.  0724 

13 

.0485 

0 

0 

.  1433 

1139 

.  1193 

14 

.0485 

0 

0 

.  3260 

2622 

.2780 

15 

0 

0 

0 

0 

0 

0 

16 

0 

0 

0 

0 

0 

0 

17 

0 

0 

0 

0 

0 

0 

18 

0 

0 

0 

0 

0 

0 

19 

0 

0 

0 

0 

0 

0 

20 

0 

0 

0 

0 

0 

0 

21 

0 

0 

0 

0 

0 

0 

22 

0 

0 

0 

0 

0 

0 

23 

0 

0 

0 

0 

0 

0 

24 

0 

0 

0 

0 

0 

0 

25 

0 

0 

0 

0 

0 

0 

26 

0 

0 

0 

0 

0 

0 

27 

0 

0 

0 

0 

0 

0 

28 

0 

0 

0 

0 

0 

0 

29 

0 

0 

0 

0 

0 

0 

30 

0 

0 

0 

0 

0 

0 
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TABLE  21 

AVIATION  ATTRITION  RATES  FOR  LTCOLS 


LOS 


AGG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

C 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  1457 

0 

0 

.2694 

.2345 

.2438 

.  1457 

.0596 

.0952 

.2517 

.2305 

.2343 

.  1457 

.0861 

0 

.2370 

.0251 

.0251 

.  1457 

.  1948 

0 

.0084 

.0106 

.  0096 

.  1457 

.2444 

.0055 

.0102 

.0131 

.0113 

.  1457 

.2669 

.  1158 

.  1221 

.  1260 

.  1252 

.  1457 

.2726 

.0928 

.  1005 

.  1050 

.  1033 

.  1457 

.2788 

.  1670 

.  1697 

.  1729 

.  1729 

.  1457 

.2667 

.2098 

.2118 

.2135 

.2135 

.  1457 

.2352 

.2012 

.2087 

.2096 

.2096 

.1457 

.2019 

.2568 

.2622 

.2602 

.2602 

.  1457 

.  1268 

.2877 

.  3038 

.2973 

.2973 

.  1457 

.0384 

.  3288 

.3368 

.  3241 

.  3260 

.  1457 

.0003 

.3415 

.  3303 

.3123 

.3161 

.  1457 

.0199 

.2963 

.  3078 

.2864 

.2911 

.  1457 

.0244 

.5385 

.5251 

.4940 

.  5066 

TABLE  22 

COMBAT  SUPPORT  ATTRITION  RATES  FOR  1ST  LTS 
CODE  07  ENGINEERS 


OS 

A GG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

2 

.2267 

.3593 

.2111 

.  1706 

.  2035 

.17  7  6 

3 

.2267 

.3989 

.2465 

.2551 

.  2817 

.2609 

4 

.2267 

.  3957 

.  1970 

.2057 

.2396 

.2112 

5 

.2267 

.3812 

.2266 

.2257 

.2540 

.2322 

6 

.2267 

.  3534 

.  3030 

.3112 

.3193 

.3193 

7 

.2267 

.3007 

.  3059 

.  3304 

.  3245 

.  3245 

8 

.2267 

.  1799 

.4444 

.4945 

.4282 

.4721 

9 

.2267 

.  1958 

.  1935 

.2290 

.2224 

.2224 

10 

.  2267 

.2064 

.  1765 

.2128 

.2115 

.2115 

11 

.2267 

.  1995 

.  1250 

.  1646 

.  1712 

.  1712 

12 

.2267 

.  1755 

.0769 

.1196 

.  1298 

.  1298 

13 

.2267 

.  1506 

0 

.0775 

.  0901 

.0901 

14 

.2267 

.  1755 

0 

.0700 

.0873 

.0821 

15 

.2267 

.  1661 

0 

.0702 

.  0861 

.  0827 

16 

.2267 

.  1263 

0 

.0891 

.0959 

.0959 

17 

.2267 

.  1194 

0 

.0910 

.0987 

.0987 

18 

.2267 

.  1506 

0 

.0760 

.0887 

.0887 

19 

.2267 

.  1263 

0 

.0891 

.0959 

.  0959 

20 

.2267 

.  1047 

.  1429 

.  1972 

.  1773 

.  1773 

21 

.2267 

.  1122 

.  1333 

.2050 

.  1852 

.  1852 

22 

.2267 

.0967 

0 

.  1151 

.  1114 

.  1114 

23 

.2267 

.0703 

0 

.  1286 

.1161 

.  1161 

24 

.2267 

.0196 

0 

.2507 

.  1889 

.2155 

25 

.2267 

.0001 

0 

.3260 

.2203 

.2793 

26 

.2267 

.0001 

0 

.  3260 

.2203 

.2793 

27 

0 

0 

0 

0 

0 

0 

28 

0 

0 

0 

0 

0 

0 

29 

0 

0 

0 

0 

0 

0 

30 

0 

0 

0 

0 

0 

0 
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TABLE  23 


.OS 


COMBAT  SUPPORT  ATTRITION  RATES  FOR  1ST  LTS 
CODE  13  COMMUNICATIONS 


GG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.2267 

.3601 

.2418 

.  1977 

.2268 

.205 

.2267 

.4080 

.2910 

.2846 

.3078 

.2901 

.2267 

.4046 

.2687 

.2664 

.2921 

.2719 

.2267 

.3910 

.3399 

.3417 

.  3511 

.3484 

.2267 

.  3563 

.2558 

.2673 

.2841 

.2757 

.2267 

.2963 

.4198 

.4319 

.4048 

.4186 

.2267 

.2879 

.2432 

.2589 

.2645 

.2645 

.2267 

.  3257 

.2281 

.2335 

.2507 

.2433 

.2267 

.  3330 

.  1120 

.  1132 

.  1489 

.  1203 

.2267 

.2940 

.0760 

.  1015 

.  1326 

.  1100 

.2267 

.2218 

.0513 

.0920 

.  1135 

.  1033 

.2267 

.  1122 

.  1333 

.  1922 

.  1753 

.  1753 

.2267 

.0403 

0 

.  1648 

.  1348 

.  1379 

.2267 

.0001 

0 

.  3629 

.2467 

.3147 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.2267 

.0001 

0 

.  3260 

.2203 

.2793 

.2267 

.0001 

0 

.  3260 

.2203 

.2793 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

TABLE  24 

COMBAT  SUPPORT  ATTRITION  RATES  FOR  1ST  LTS 
CODE  20  MOTOR  TRANSPORT 


AGG  ORIG  AGG  TRANS 


0 

0 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 

.2267 


0 

0 

.  1661 
.2508 
.2724 
.2784 
.2189 
.0796 
.0299 
.0607 
.0403 
.  1047 
.  1047 
.  1506 
.  1449 
.  1841 
.  1390 
.0884 
.0703 
.0796 
.0507 
.0703 
.0403 
.0403 
.0103 
.0103 
.0103 
.0029 
.  0063 
.0063 
.0103 


0 

0 

0 

3137 

2857 

1791 

2632 

1818 

1 

6667 

0 

1429 

0 

0 

0 

0 

0 

0 

0 

0 

2500 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


TSCA 

0 

0 

0878 

4064 

2693 

2207 

3449 

3036 

8007 

6184 

1866 

1935 

1094 

0810 

0897 

0651 

0816 

1146 

1343 

1293 

3191 

1343 

2273 

1866 

2748 

2687 

2687 

2632 

4033 

4033 

2394 


0 

0 

1014 

3747 

2699 

2316 

3191 

2517 

6388 

4918 

151 

1745 

1084 

0931 

0996 

0841 

0917 

1093 

1204 

1188 

2535 

1204 

1811 

1510 

2004 

1960 

1960 

1848 

2631 

2631 

1749 


LTJS 

0 

0 

1014 

3899 

2699 

2316 

3267 

2744 

7683 

5845 

1582 

1745 

1084 

0931 

0996 

0764 

0917 

1093 

1204 

1188 

2860 

1204 

1966 

1582 

2364 

2305 

2305 

2227 

3492 

3492 

2028 


TABLE  25 


COMBAT  SUPPORT  ATTRITION  RATES  FOR  LTCOLS 
CODE  07  ENGINEERS 


AGG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  . 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  I486 

.0253 

0 

.3629 

.2847 

.2847 

0 

0 

0 

0 

0 

0 

.  1486 

.  1767 

0 

.0967 

.  1090 

.  1090 

.  I486 

.2900 

0 

.0392 

.0675 

.0531 

.  1486 

.3117 

0 

.0323 

.0620 

.0436 

.  1486 

.3184 

.2500 

.2413 

.2541 

.2541 

.  1486 

.3158 

.0645 

.0888 

.  1197 

.  1060 

.  1486 

.3088 

.  1404 

.  1387 

.  1643 

.  1601 

.  1486 

.2860 

.0909 

.  1171 

.  1419 

.  1397 

.  1486 

.2586 

.2424 

.2597 

.2595 

.2595 

.  1486 

.2219 

.2609 

.2922 

.2797 

.2797 

.  1486 

.  1244 

0 

.  1433 

.  1400 

.  1400 

.  1486 

.1132 

.  7500 

.  6913 

.5861 

.62  72 

.  1486 

.0420 

.  6667 

.  5462 

.4389 

.4577 

.  1486 

.0420 

0 

.  2954 

.2394 

.2394 

.  1486 

.  1010 

.5714 

.4929 

.4150 

.4224 

TABLE  26 


COMBAT  SUPPORT  ATTRITION  RATES  FOR  LTCOLS 
CODE  13  COMMUNICATIONS 


LOS 

AGG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

0 

5 

0 

0 

0 

0 

0 

0 

6 

0 

0 

0 

0 

0 

0 

7 

0 

0 

0 

0 

0 

0 

8 

0 

0 

0 

0 

0 

0 

9 

0 

0 

0 

0 

0 

0 

10 

0 

0 

0 

0 

0 

0 

11 

0 

0 

0 

0 

0 

0 

12 

0 

0 

0 

0 

0 

0 

13 

0 

0 

0 

0 

0 

0 

14 

0 

0 

0 

0 

0 

0 

15 

0 

0 

0 

0 

0 

0 

16 

.  1486 

.01 

0 

.4033 

.  3067 

.3067 

17 

.  1486 

.2018 

0 

.  1091 

.  1233 

.  1233 

18 

.  1486 

.  2817 

0 

.0445 

.0723 

.0598 

19 

.  1486 

.3221 

0 

.02955 

.0596 

.  0396 

20 

.  1486 

.3378 

.  1463 

.  1762 

.2014 

.  1959 

21 

.  1486 

.3330 

.  1299 

.  1713 

.  1964 

.  1913 

22 

.  1486 

.  3330 

.  1558 

.  1910 

.2135 

.2119 

23 

.  1486 

.  3026 

.  1509 

.  1491 

.  1726 

.  1719 

24 

.  1486 

.2817 

.2857 

.2577 

.2617 

.2617 

25 

.  1486 

.2524 

.  3226 

.2887 

.2823 

.2823 

26 

.  1486 

.  1900 

.2353 

.2658 

.2522 

.2522 

27 

.  1486 

.  1444 

.3636 

.3866 

.  3401 

.  3401 

28 

.  1486 

.  1132 

0 

.  1775 

.  1656 

.  1656 

29 

.  1486 

.  1132 

.2500 

.  3483 

.  3019 

.  3019 

30 

.  1486 

.2586 

.  3636 

.3884 

.  3652 

.  3652 

TABLE 

27 

COMBAT  SUPPORT  ATTRITION  RATES  FOR 

LTCOLS 

CODE  20  MOTOR 

TRANSPORT 

.OS 

AGG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LT  JS 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

0 

5 

0 

0 

0 

0 

0 

0 

6 

0 

0 

0 

0 

0 

0 

7 

0 

0 

0 

0 

0 

0 

8 

0 

0 

0 

0 

0 

0 

9 

0 

0 

0 

0 

0 

0 

10 

0 

0 

0 

0 

0 

0 

11 

0 

0 

0 

0 

0 

0 

12 

0 

0 

0 

0 

0 

0 

13 

0 

0 

0 

0 

0 

0 

14 

0 

0 

0 

0 

0 

0 

15 

0 

0 

0 

0 

0 

0 

16 

.  1486 

.0100 

0 

.4033  . 

3067 

.  3067 

17 

.  1486 

.0100 

0 

.4033  . 

3067 

.  3067 

18 

.  1486 

.0253 

0 

.3260  . 

2562 

.2562 

19 

.  1486 

.0558 

0 

.2113  . 

1794 

.  1794 

20 

.  1486 

.0878 

3333 

.3641  . 

3078 

.3078 

21 

.  1486 

.0736 

4000 

.4277  . 

3542 

.3542 

22 

.  1486 

.0736 

4000 

.4277  . 

3542 

.  3542 

23 

.  1486 

.0100 

0 

.4033  . 

3067 

.  3067 

24 

0 

0 

0 

0 

0 

0 

25 

0 

0 

0 

0 

0 

0 

26 

0 

0 

0 

0 

0 

0 

27 

0 

0 

0 

0 

0 

0 

28 

0 

0 

0 

0 

0 

0 

29 

.  1486 

.0253 

0 

.3260  . 

2562 

.2562 

30 

.  1486 

.0878 

6667 

.6067  . 

5057 

.5333 

TABLE  28 

GROUND  COMBAI 

ATTRITION  RATES  FOR 

1ST 

_,TS 

CODE  03 

INFANTRY 

OS 

AGG  ORIG  AG 

G  TRANS 

MLE 

TSCA 

JS 

LTJS 

0 

0 

0 

0 

0 

0 

0 

1 

.  1991 

0 

0 

.4033 

.2713 

.3177 

2 

.1991 

3347 

.  1788 

.  1382 

.  1514 

.  1438 

3 

.1991 

3814 

.  1935 

.  1899 

.2034 

.  1944 

4 

.1991 

3815 

.2165 

.2148 

.2269 

.2195 

5 

.1991 

3709 

.  1746 

.  1748 

.  1885 

.  1796 

6 

.1991 

3434 

.2604 

.2581 

.2646 

.2646 

7 

.1991 

2973 

.2194 

.2250 

.2305 

.2305 

8 

.1991 

2283 

.  1463 

.  1497 

.  1554 

.  1554 

9 

.1991 

2626 

.  1511 

.  1548 

.  1626 

.  1626 

10 

.1991 

2622 

.  1516 

.  1611 

.  1685 

.  1685 

11 

.1991 

2241 

.  0808 

.0920 

.  1006 

.  1002 

12 

.1991 

1013 

.  1039 

.  1240 

.  1222 

.  1222 

13 

.1991 

0186 

.0541 

.0928 

.  0851 

.0850 

14 

.1991 

0626 

0 

.  1179 

.0887 

.0887 

15 

.1991 

0626 

.  1667 

.2294 

.  1811 

.  1889 

16 

.1991 

2026 

0 

.  1913 

.  1372 

.  1463 

17 

0 

0 

0 

0 

0 

0 

18 

0 

0 

0 

0 

0 

0 

19 

0 

0 

0 

0 

0 

0 

20 

0 

0 

0 

0 

0 

0 

21 

0 

0 

0 

0 

0 

0 

22 

0 

0 

0 

0 

0 

0 

23 

0 

0 

0 

0 

0 

0 

24 

0 

0 

0 

0 

0 

0 

25 

0 

0 

0 

0 

0 

0 

26 

0 

0 

0 

0 

0 

0 

27 

0 

0 

0 

0 

0 

0 

28 

0 

0 

0 

0 

0 

0 

29 

0 

0 

0 

0 

0 

0 

30 

0 

0 

0 

0 

0 

0 
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TABLE  29 

GROUND  COMBAT  ATTRITION  RATES  FOR  1ST  LTS 
CODE  05  ARTILLERY 


AGG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  1991 

.2888 

.  1170 

.  1050 

.  1168 

.  1115 

.  1991 

.3408 

.  2052 

.2096 

.2193 

.2160 

.  1991 

.3423 

.2462 

.2480 

.2552 

.2547 

.  1991 

.3197 

.  1948 

.  1992 

.2080 

.2063 

.  1991 

.2724 

.3016 

.3049 

.3023 

.  3023 

.  1991 

.  1671 

.3333 

.  3329 

.3184 

.3184 

.  1991 

.  1294 

.  1263 

.  1323 

.  1321 

.  1321 

.  1991 

.  1390 

.2353 

.2452 

.2360 

.  2  "  .0 

.  1991 

.  1503 

.0901 

.0991 

.  1029 

.  1029 

.  1991 

.0888 

.0571 

.0783 

.0791 

.0791 

.  1991 

.0342 

.0454 

.0851 

.0803 

.0803 

.  1991 

.  0039 

.  1000 

.  1609 

.  1341 

.  1341 

.  1991 

.2538 

0 

.  1758 

.  1216 

.  1305 

.  1991 

0 

0 

.  3260 

.2189 

.2523 

.  1991 

0 

0 

.3260 

.2189 

.2523 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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TABLE  30 

GROUND  COMBAT  ATTRITION  RATES  FOR  1ST  LTS 
CODE  10  TANKS 


;g  orig 

AG G  TRANS 

MLE 

TSCA 

JS 

LTJS 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1991 

.  1573 

.2906 

.2368 

.2300 

.2300 

1991 

.2506 

.2490 

.2505 

.  2506 

.2506 

1991 

.2515 

.2390 

.2435 

.2441 

.2441 

1991 

.2488 

.  1633 

.  1697 

.  1756 

.  1756 

1991 

.  1987 

.  3230 

.  3216 

.3112 

.3112 

1991 

.  1030 

.3333 

’.3449 

.  3223 

.  3241 

1991 

.0186 

.  1081 

.  1377 

.  1246 

.  1246 

1991 

.0364 

.  1778 

.2113 

.  1931 

.  1931 

1991 

.0651 

.  1034 

.  1167 

.  1121 

.  1121 

1991 

.0207 

.0526 

.0855 

.  0789 

.  0789 

1991 

.0801 

0 

.  1498 

.  1131 

.  1149 

1991 

0 

0 

.  1927 

.  1311 

.  1431 

1991 

0 

0 

.4033 

.2713 

.3177 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1991 

0 

0 

.  3260 

.2189 

.2523 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

TABLE  31 

GROUND 

COMBAT  ATTRITION  RATES  FOR 

LTCOLS 

CODE  03 

INFANTRY 

AGG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LT  JS 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  1278 

.  1361 

0 

.  3002 

.2664 

2749 

.  1278 

.0249 

0 

.  1000 

.0962 

0962 

.  1278 

.2097 

0 

.0173 

.0210 

0193 

.  1278 

.3101 

0 

.0062 

.  0098 

0070 

.  1278 

.  3341 

.0095 

.0139 

.0188 

0148 

.  1278 

.3431 

.1224 

.  1273 

.  1342 

1299 

.  1278 

.  3395 

.0752 

.0804 

.0877 

0826 

.  1278 

.  3328 

.  1253 

.  1285 

.  1351 

1313 

.  1278 

.3155 

.  1365 

.  1413 

.  1471 

1445 

.  1278 

.  3001 

.2113 

.2171 

.2202 

2202 

.  1278 

.2658 

.2090 

.2163 

.2182 

2182 

.  1278 

.2236 

.2464 

.2536 

.2524 

2524 

.  1278 

.  1703 

.2889 

.2964 

.2911 

2911 

.  1278 

.0984 

.2308 

.2343 

.2281 

2281 

.  1278 

.  0744 

.2326 

.2111 

.2047 

2047 

.  1278 

.  1550 

.4500 

.4551 

.4421 

4460 

TABLE  32 


GROUND  COMBAT  ATTRITION  RATES  FOR  LTCOLS 
CODE  05  ARTILLERY 


G  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1278 

.0215 

0 

.  1994 

.  1807 

.  1818 

1278 

.  1255 

0 

.  0314 

.0340 

.0340 

1278 

.2219 

0 

.0141 

.0177 

.0158 

1278 

.2496 

.0347 

.  0477 

.0528 

.0504 

1278 

.2463 

.2500 

.2530 

.2527 

.2527 

1278 

.2357 

.0654 

.0747 

.0795 

.0782 

1278 

.  2201 

.  1045 

.0914 

.0956 

.0956 

1278 

.  1826 

.  1414 

.  1622 

.  1630 

.  1630 

1278 

.  1533 

.  3544 

.  3577 

.  3487 

.  3488 

1278 

.0984 

.  1923 

.2046 

.  1999 

.  1999 

1278 

.0597 

.5263 

.5180 

.4962 

.5050 

1278 

.0329 

.6667 

.5840 

.5432 

.5610 

1278 

.0215 

.4000 

.4096 

.  3774 

.  3876 

1278 

.  1361 

.  8000 

.  6480 

.  5981 

.  6210 

1278 

.0133 

.  7273 

.6592 

.6192 

.  6385 

TABLE  33 


GROUND  COMBAT  ATTRITION  RATES  FOR  LTCOLS 


CODE  10 

TANKS 

AGG  ORIG 

AGG  TRANS 

MLE 

TSCA 

JS 

LTJS 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  1278 

0 

0 

.3260 

.2811 

.2943 

.  1278 

.0037 

0 

.  1354 

.  1239 

.  1239 

.  1278 

.0129 

0 

.0695 

.  0665 

.0665 

.  1278 

.0716 

0 

.0427 

.0437 

.0437 

.  1278 

.0744 

.  1860 

.  1828 

.  1778 

.  1778 

.  1278 

.0628 

.0513 

.0861 

.0851 

.0851 

.  1278 

.0773 

.  1818 

.2047 

.  1988 

.  1988 

.  1278 

.0567 

.0541 

.0883 

.0870 

.0870 

.  1278 

.0505 

.2857 

.2982 

.2857 

.2860 

.  1278 

.0129 

.3478 

.4256 

.4026 

.4096 

.  1278 

.0002 

.  1333 

.2135 

.  1978 

.  1980 

.  1278 

.0037 

.  3077 

.3171 

.2933 

.2983 

.  1278 

.0329 

.2222 

.2630 

.2385 

.2428 

.  1278 

.0699 

.2857 

.  3350 

.  3027 

.3113 

.  1278 

.0001 

.  3750 

.  3841 

.  3590 

.  3659 

V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

This  study  investigated  the  performance  of  various 
attrition  rate  estimation  schemes,  with  particular  emphasis 
on  small  cell  performance.  The  hope  was  to  identify  a 
specific  method,  the  limited  translation  James-Stein  tech¬ 
nique,  as  the  best  predictor  of  cell  rates.  This  has  not 
been  done.  In  fact,  the  relative  performance  of  four  schemes 
(MLE,  TSCA,  James-Stein,  limited  translation  James-Stein) 
was  so  varied  that  no  clear  winner  or  pattern  of  performance 
emerged . 

Several  points  are  of  interest.  First,  even  though  the 
performance  of  the  four  techniques  listed  above  varied,  they 
were  all  uniformly  better  than  either  original  or  trans¬ 
formed  scale  aggregate  estimation.  Thus  four  excellent 
candidates  for  replacing  the  present  attrition  rate 
estimation  methods  are  available  for  further  testing. 

Second,  the  models  seem  very  sensitive  to  small  changes 
in  parameters.  When  investigating  the  behavior  of  the 
Freeman-Tukey  transform,  and  James-Stein  and  limited  trans¬ 
lation  James-Stein  estimation,  several  choices  had  to  be 
made  concerning  the  methods  used  to  calculate  various  quan¬ 
tities,  e.g.,  the  grand  mean,  maximum  likelihood  estimators, 
and  the  inventories.  Large  changes  in  FOM  values  and  cell 
attrition  rates  were  observed  as  different  methods  were 
tried.  This  lack  of  robustness  is  troublesome,  and  indicates 
the  models  are  not  ready  for  implementation.  We  feel  a 
better  method  of  aggregating  the  OF/LOS/grade  cells  would  do 
much  to  relieve  this  problem. 

Third,  in  the  small  cells,  the  Freeman-Tukey  transforma¬ 
tion  fails  to  normalize  the  cell  means  or  stabilize  the 
variance.  This  failure  is  an  inherent  aspect  of  dealing  with 


small  cells.  In  general,  the  research  effort  must  be  alert 
to  finding  alternative  ways  to  manage  the  small  cell 
problem. 


B.  RECOMMENDATIONS 

At  present,  no  method  investigated  here  is  recommended 
for  implementation.  Further  study  in  the  following  areas  is 
needed . 

1.  Aggregation.  The  work  of  Amin  Elseramegy  [Ref.  2] 
needs  to  be  carried  forward  to  identify  statistically 
well  behaved  aggregates. 

2.  TSCA.  The  performance  of  the  TSCA  scheme  should  be 
investigated  in  light  of  its  surprising  performance. 

3.  Yearly  Update.  The  problem  of  when  and  how  to  update 
the  estimation  data  base  must  be  solved. 

4.  Use  of  Different  Estimators.  Investigate  the  use  of 
several  different  estimators  (such  as  MLE ,  TSCA, 
James-Stein,  and  limited  translation  James  -  Stein) , 
either  in  combination  to  yield  average  cell  estimates, 
or  separately  to  estimate  the  rates  of  different  cells 
based  on  OF,  military  occupation  specialty  (MOS)  as  a 
subset  of  OF,  LOS,  grade,  or  desired  estimation  year. 
This  would  require  tne  aggregation  problem  solved. 

5.  Robust  Parametric  Empirical  Bayes.  Investigate  the  use 
of  the  robust  parametric  empirical  Bayes  model  to 
determine  global  and  small  cell  estimation  efficacy. 


APPENDIX  A 

ESTIMATION  ALGORITHMS 


1.  JAMES -STEIN  ALGORITHM 

As  used  by  Major  Tucker  [Ref.  1],  this  algorithm  calcu¬ 
lates  the  James-Stein  estimator  of  attrition  rates.  The 
following  formulas  use  a  double  indexing  system,  (i,j),  to 
identify  cells  in  an  aggregate  and  use  t  =  l,...,n  to  iden¬ 
tify  time  periods.  This  usage  is  different  from  that  of  the 
main  text  which  used  the  single  subscript  i  =  l,...,k  to 

identify  cells.  This  was  convenient  for  purposes  of  ANOVA. 
The  current  double  indexing  is  needed  to  access  the  real 
data  in  its  natural  form. 

Notation: 

I  =  number  of  LOS  cells  in  the  chosen  aggregate 
J  =  number  of  OF  cells  in  the  chosen  aggregate 
INVjj(t)  =  inventory  with  LOS  =  i  and  OF  =  j  at 
beginning  of  year  t,  t=l,...,T 
y^j(t)  =  number  of  attritions  in  cell  (i,j)  during 
year  t 

nij(t)  =  maximum  { y i j ( t ) , ■ 5 [INVi j ( t ) + INVi j ( t + 1 ) ] } 

D  -  incidence  matrix  which  identifies  cells  with 

sampling  zeros  (average  inventory  zero  for  all 
estimation  years) 

=  0  if  cell  is  a  structural  zero 
=  1  if  cell  is  not  a  structural  zero 


STEP  1:  Use  a  variance  stabilizing  transformation 
( Freeman-Tukey ) . 


xij(t)  =  0 . 5[ni i (t)+0. 5]^{sin' 1[2(y- i (t)/ (ni i (t )+l)-l]  (A.l) 


U 


sin'1[2(yii (t)  +  l)/ (nii (t)  +  l)-l]  )  . 


ij 


ij 


STEP  2:  Form  the  cell  means  and  the  grand  mean. 


Xij  =  (1/T)pij(t),  V 

t 


(A. 2) 


X  =  (l/K)£pijDij,  V  i,J, 

'  j 


(A. 3) 


K  =  IZDij’  V  i,j.  (A. 4) 

'  j 

STEP  3:  Form  the  sum  of  squares  error  (SSE),  and  the  sum 
of  squares  between  (SSB),  using  the  total  sum  of  squares 
( SST )  . 


SST  =  IXXUijCtO-xPDij,  V  i,j,t 

‘  J  t 


(A. 5) 


SSE  =  IH[xij(t)-xij]2Dij,  V  i,j,t  (A. 6) 

»  J  t 

SSB  =  SST- SSE  (A. 7) 


STEP  4:  Compute  the  set  of  James-Stein  estimators  in  the 
transformed  scale. 

Cj  =  1  -  (K-3)SSE/ (K(T-1)+2)SSB  (A. 8) 


X  ♦  Cj(xirx) 

Dij  =  l 

Pj(i,J)  = 

if 

undefined 

D,-  i  -  0 

(A. 9) 


STEP  5:  Invert  the  transform  to  produce  the  attrition 
rates  r^j .  This  corrects  a  typographical  error  that  appeared 
in  this  step  in  reference  1. 


nij  =  (l/T)Xnij(t) 


vtj  =  Pj(i, j )/ (nij+ .5)^ 


vij  <  - 7t /  2 


(A. 10) 


(A. 11) 


(A. 12) 


. 5 [1+sin  ytj] 


-7l/2<Vi  j<7l/2" 
vij  >n/2 


2.  LIMITED  TRANSLATION  JAMES- STEIN  ALGORITHM 

This  algorithm  calculates  the  limited  translation 
James-Stein  estimator  of  attrition  rates.  Because  the  only 
difference  between  this  algorithm  and  the  James-Stein  algo¬ 
rithm  of  Appendix  A,  Section  1  is  in  STEP  4,  only  that 
changed  step  is  presented.  All  other  steps  remain  the  same. 


STEP 


Compute  the  set  of  limited  translation 


James-Stein . estimators  in  the  transformed  scale.  Choose  d 
from  the  range  [0,oo].  Values  between  .2  and  1.0  seem  to 
give  the  best  results. 


1  -  (K-3)SSE/(K(T-1)+2)SSB 


(K(T-1)+2)SSB/ (K-3)SSE  -  1 


(A. 13) 


(A. 14) 


LI  =  (X±j  -X)2/  (A+l) 

p(u)  =  minimum  (l.d/u'2) 


CLJ  =  1  -  p(u)[(K-3)SSE/(K(T-l)+2)SSB] 


x  ♦  cLJ(xirx) 


Dij: 


^LJ ( i ’ J )  " 


if 


undefined 


D,  ,= 


APPENDIX  B 

UTILIZATION  OF  THE  FREEMAN- TUKEY  ARCSINE  TRANSFORMATION 


1.  GENERAL  I 

Limited  translation  James-Stein  estimation  as  discussed 
in  the  Efron  and  Morris  articles  [Ref.  5,6]  makes  the 
following  assumptions: 

(1)  The  distribution  of  the  number  of  losses  is  normally 
distributed . 

(2)  The  variances  of  the  normally  distributed  number  of 
losses  are  equal. 

Since  the  raw  losses  are  not  normally  distributed,  the 
Freeman-Tukey  arcsine  transformation 

x  =  0 . 5  [n+0 . 5]^{  sin"  [2y/ (n+ 1)  -  1]  (B.l) 

+  sin" ^ [2 (y+ 1 ) / (n+ 1 ) -  1] } 


was  used,  given  the  central  inventory  n,  to  transform  the 
raw  losses  y  into  transformed  losses  x. 

2.  NORMALITY  OF  TRANSFORMED  LOSSES 

Tucker  [Ref.  1]  demonstrated  that  the  distribution  of  y 
given  n  is  Binomial  (p),  where  the  parameter  p  is  the  prob¬ 
ability  of  an  individual  loss.  The  probability  mass  function 
of  this  discrete  distribution  is 

P( Y  =  y }  =  [n!  /  (y!  (n-y)  !  )]  py(l-p)n'y.  (B.2) 

Since  there  is  a  unique  mapping  from  Y  to  X,  it  is  true  that 


P{X=x}  =  P{ Y=y } 


(B .  3  ) 


To  check  normality,  a  central  inventory  n  was  chosen. 
Possible  losses  y,  from  0  to  n,  were  identified.  Then  the 
above  probability  mass  function  and  transformed  losses  x 
were  calculated.  This  data  was  then  displayed  on  an  x  versus 
P{X=x}  graph.  Figures  B.l,  B.2,  and  B.3  are  graphs  for  n 
equals  5,  8,  and  10.  Probabilities  of  loss  are  for  p  equals 
.05,  .1,  .15,  and  .2.  The  ranges  of  greatest  concern  were 
those  of  low  central  inventory  n  and  low  probability  of  loss 
P- 

By  inspection,  the  higher  the  values  of  n  and  p,  the 
better  is  the  normal  approximation.  The  Freeman-Tukey 
transformation  is  unreliable  at  low  values  of  n  and  p. 

3.  VARIANCE  STABILITY  OF  TRANSFORMED  LOSSES 

The  mean  and  variance  of  the  Binomial  distribution  is 
well  known.  Again  using  the  fact  that  if  x  is  a  function  of 
y  and  n,  then 

P{X=x}  =  P{ Y=y } ,  (B . 5 ) 

and  the  variance  of  x  can  be  calculated,  given  values  of  n, 
y  and  x.  The  values  of  the  variance  of  x  were  calculated, 
and  are  graphed  in  Figure  B.4. 

By  inspection,  once  n  equals  7,  the  variance  stability 
at  values  of  p  less  than  .2  is  poor.  Again,  the 
Freeman-Tukey  transformation  is  unreliable  at  low  values  of 
n  and  p . 


APPENDIX  C 

LIMITED  TRANSLATION  FACTOR  DERIVATION 


This  is  a  derivation  of  the  form  of  the  limited  transla 
tion  factor 

p ( u )  =  minimum  (l,d/u^)  (C.l) 

Given:  From  Efron  and  Morris  [Ref.  5], 

x+M  x<-C  (C.2) 

SAiM(x)  =  Ax/(A+1)  if  xs [ -C , C] 

x-M  x>C 

where,  if  Xjj  and  X  are  those  used  in  Appendix  A, 

C  =  M(A+1)  ( C . 3 ) 


d  =  M(A+1)* 


(C.4) 


u  =  (X± j -X)2  /  ( A+ 1 ) 


(C.5) 


5a  m(x)  =  [l-p(u)/ (A+ I)]x. 


(C.6) 


Derivation:  From  equation  C.2  we  have, 


x+M  x<-C  (C . 7 ) 

SAjM(x)  =  Ax/ (A+l)  if  xe [-C , C] 

x-M  x>C 

By  substitution,  and  multiplication  by  x/x, 

(x+M)x/x  [u (A+ 1 ) ]^< -M(A+ 1 )  (C.8) 

(l-p(u)/ (A+l)  )x  =  Ax/(A+1)  if  [u(A+l)]'se  [-M(A+1)  ,M(A+1)] 

(x-M)x/x  [u(A+ 1)]'S>M(A+1) 

By  multiplication  by  1/x,  and  cancellation, 

(x+M)x/x  u^< -M(A+ 1 )^  (C.9) 

1  - p ( u ) / (A+l)  =  Ax/ ( A+ 1 )  if  u^c  [-M(A+1)'^,M(A+1)'S] 

(x-M)x/x  u's>M(A+l)'s 

By  subtracting  1,  and  multiplying  by  -1, 

1- ( x+M) / x  u^< -M(A+ 1 (C.10) 

p(u)/(A+l)  =  1-A/(A+1)  if  u^e  [-M(A+1)'S,M(A+1)'^] 
l-(x-M)/x  u's>M(A+l)'s 


By  multiplying  by  (A+l), 

( A+ 1 ) [ 1- (x+M) / x] 
p ( u )  =  (A+l) [1-A/ (A+l)] 

(A+l) [1- (x-M)/x] 


-M( (A+l)/u)^ 
1 

M ( (A+l)/u)* 


and  substitution, 

u%<-d  (C.ll) 

if  u^e  [-d,d] 

At  . 

u  ->d 


u %<-d  (C . 12 ) 

if  u^e[-d,d] 
u^>d 


If  u^<-d,  then  0<-d/u'£<l,  and  if  u^>d,  then  0<d/u's<l. 
Therefore,  the  last  equation  may  be  written  as 


p(u)  =  minimum  (l,d/u'£). 

It  is  worth  noting  that  if  xg[-C,C]  and 

SAjM(x)  =  Ax/(A+1). 

Also,  if  x^[-C,C]  and  0<p(u)<l,  then 

Ax/(A+1)  <  5a?m(x)  <  x. 

Thus,  if  x  is  outside  the  interval  [~C,C] 
translated  less  than  any  x's  inside  the 
This  confirms  the  graphical  representat 
translation  shown  in  Figure  3.1. 


(C.14) 

p(u)=l,  then 

(C.15) 

(C.16) 

,  then  this  x  is 
interval  [  - C , C] . 
ion  of  limited 


APPENDIX  D 
DATA  MANIPULATION 

1 .  GENERAL 

This  appendix  documents  the  steps  and  programs  used  to 
take  the  raw  data  from  magnetic  tape  to  the  output  of 
figures  of  merit  and  attrition  rates.  The  program  languages 
include  JCL,WATFIV,  and  APL.  The  author  is  familiar  with  the 
system  at  the  Naval  Postgraduate  School,  and  the  instruc¬ 
tions  are  therefore  specific  for  that  system.  However,  minor 
changes  should  be  all  that  are  necessary  to  implement  these 
procedures  elsewhere.  To  achieve  good  results,  the 
procedures  should  be  followed  in  the  order  presented. 

2.  CONVERSION  OF  RAW  DATA  FROM  TAPE  TO  APL 

The  original  data  is  on  a  magnetic  tape  named  COUNTS, 
prepared  by  NPRDC .  The  tape  is  held  by  either  Professor 
R.R.  Read,  or  by  the  Computer  Center  in  Professor  Read's 
name.  Ensure  the  tape  is  properly  logged  into  the  Computer 
Center,  and  submit  the  JCL  program  IEBGENER  in  Figure  D.I  to 
put  the  tape  on  mass  storage.  The  data  set  in  mass  storage 
should  be  named  MSS . SXXXX . COUNTS ,  where  XXXX  is  the  user  ID 
number  of  the  operator. 

Once  in  mass  storage,  submit  the  JCL  program  MSSCOUN  in 
Figure  D.2  to  move  the  data  from  mass  storage  to  the  MVS004 
disk.  Data  on  this  disk  is  accessible  from  CMS. 

Use  the  system  exec  GETMVS  to  move  the  data  from  the 
MVS004  disk  to  your  disk.  Simply  enter  'GETMVS'  on  the 
computer  and  follow  the  directions.  The  identification 
requested  by  the  prompts  of  GETMVS  will  be  'SXXXX  COUNTS'. 
Since  the  data  set  is  large  (16,093  lines  of  53  columns 
each),  it  is  advisable  to  get  additional  woikspace  by  either 
applying  for  a  B  disk  of  at  least  8  cylinders,  or  by  getting 


//  EXEC  PGM=IEBGENER 
//SYSPRINT  DD  SYS0UT=A 


//SYSIN 
/ / SYSUT1 

// 

// 

/ / SYSUT2 

// 

// 

// 


DD 

DD 


DD 


DUMMY 

UNIT=3400-5 ,VOL=SER= COUNTS ,DISP= (OLD, PASS 
DCB= (RECFM=FB ,LRECL=53 , BLKSIZE=21200 , DEN= 
LABEL=7l.SL, , IN ) ,DSN= COUNTS 
UNIT  =  3330V, M$VGP  =  PUB4B , DI SP  = ( NEW . CATLG ) , 
DCB= (RECFM=FB , LRECL=53 , BLKSIZE= 1298 5 ) , 
DSN=MSS . S1662 . COUNTS 


Figure  D.l  JCL  Program 


1 


Figure  D.2  JCL  Program  MSSCOUN 

a  temporary  C  disk  of  at  least  8  cylinders.  Be  aware  that  a 
C  disk  disappears  once  logged  off,  and  all  data  on  it  is 
erased . 

The  WATFIV  program  SORT  in  Figure  D.3  should  now  be  used 
to  separate  the  data  in  COUNTS  into  seven  files,  one  file 
for  each  fiscal  year.  The  data  files  can  be  conveniently 
named  COUNXX  DATA,  where  XX  is  the  year,  e.g.,  77. 

Using  the  WATFIV  program  INV  in  Figure  D.4  ,  create  an 
array  of  inventory  indices  for  each  fiscal  year.  The  program 
should  read  in  the  data  sets  prepared  in  step  4  above,  e.g., 
COUN77  DATA.  The  output  is  read  into  a  file  that  can  be 
named  INVXX  DATA,  where  XX  is  the  fiscal  year,  e.g.,  77. 
Note  that,  for  each  fiscal  year,  a  different  DO  loop  is 
used,  since  the  number  of  data  rows  vary  form  year-to-year. 


//COUNTS  JOB  (1662,9999) , 'MAJ  J.R. ROBINSON  ,CLASS=B 
//  EXEC  PGM= IEBGENER 
//SYSPRINT  DD  SYS0UT=A 
//SYSIN  DD  DUMMY 

/ / SYSUT1  DD  DSN=MSS . S1662 . COUNTS ,DISP=SHR 
/ / SYSUT2  DD  DSN=S1662 . COUNTS ,V0L=SER=MVS004 ,UNIT=3350 , 
//  SPACE= (CYL, (2,2) , RLSE ) ,DISP= (NEW, KEEP) 


JOB 


THIS  PROGRAM  ASSIGNS  THE  DATA  FROM  COUNTS 
INTO  SEPARATE  FILES  FOR  EACH  FISCAL  YEAR. 


INV , LI , L2 , L3 , L4 , 

OF, GRADE, LOS, INV, 

8)  YEAR , OF , GRADE , 
L7,L8 

9}  YEAR , OF , GRADE , 
L7.L8 

10)  YEAR, OF, GRADE 
,L7 ,L8 

11)  YEAR, OF, GRADE 
, L7 , L8 

12)  YEAR , OF , GRADE 
,L7,L8 

13)  YEAR , OF , GRADE 
,L7 ,L8 

14)  YEAR, OF, GRADE 
,L7,L8 


INTEGER  YEAR, OF, GRADE, LOS , 
*L5 , L6 ,L7 ,L8 
DO  106  1=1,16093 
READ  (15 , 20 , END=200 )  YEAR, 
*L1,L2,L3,L4,L5 ,L6 ,L7 ,L8 
20  FORMAT  7412  ’  9l5  ) 

IF  (YEAR. EQ. 77J  WRITE  (30, 
*L0S,INV,L1,L2 .L3 ,L4.L5 ,L6 , 
IF  (YEAR.EQ.78j  WRITE  (30, 
*LOS . INV ,L1 ,L2 , L3 , L4  ,L5 .L6 , 
IF  (YEAR.EQ.79j  WRITE  [30, 

*  LOS,INV,LI,L2,L3 ,L4,L5,l8 
IF  (YEAR. EQ. 80 )  WRITE  (30, 

* , LOS , INV , LI , L2 , L3 . L4 , L5 , LS 
IF  TfEAR.EQ.8l)  WRITS  (30, 
* ,LOS,INV,LI,L2.L3 ,L4 ,L5,l£ 
IF  (YEAR. EO. 82)  WRITE  (SO, 

*  LOS . INV , LI , L2 , L3.L4.L5 , l8 
IF  (YEAR. EO. 83)  WRITE  (SO, 

*  LOS , INV , LI , L2 ,L3 , L4 , L5 , LS 
30  SormAt  (412  9lS) 

100  CONTINUE 
200  STOP 
END 

SENTRY 


Figure  D.3  WATFIV  Program  SORT 

Ensure  the  read  and  write  files  are  properly  defined  for 
each  year. 

Use  the  W7TFIV  program  LOSS  in  Figure  D.5  to  create  an 
array  of  loss  indices  for  each  fiscal  year,  similar  to  the 
inventory  indices  in  step  5  above.  Note  that  the  losses  are 
aggregated,  in  that  the  loss  data  on  the  COUNXX  DATA  file  is 
broken  into  8  different  loss  categories  (see  Reference  1). 
Such  a  breakdown  is  not  used  in  the  present  work.  See  Tucker 
[Ref.  1]  for  programs  and  procedures  if  loss  type  data  is 
desired.  Again,  the  data  can  be  read  into  files  named  LOSSXX 
DATA,  where  XX  is  the  fiscal  year.  The  arrays  LOSSXX  will  be 
significantly  smaller  than  the  INVXX  arrays. 

Finally,  use  the  APL  system  exec  CMSIO  to  move  the  data 
files  INVXX  and  LOSSXX  into  an  APL  workspace.  The  APL  arrays 


$  JOB 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


20 


100 

200 

C 

C 

C 

C 


300 

SENTRY 


INTEGER  GRADE , OF , YEAR , LOS , INDEX , IN , TOT , I YR 

THE  DO  LOOPS  CORRESPOND  TO  THE  NUMBER  OF 
RECORDS  FOR  YEARS  1977  THRU  1983. 

DO  100  1=1,2203 
DO  100  1=1,2231 
DO  100  1=1,2337 
DO  100  1=1,2351 
DO  100  1=1  2317 
DO  100  1=1,2324 
DO  100  1=1.2330 

READ  (10.20,END=200)  YEAR , OF , GRADE , LOS , INV 
•  ,L1,L2 ,L3 .L4.L5 ,L6 ,L7 ,L8 
k)RftAT  (4I2.$I5) 

IYR  =  YEAR- 76 
TOT  =  100000000*IYR 

NV.GT.O)  CALL  SUM  ( IN , TOT , OF , GRADE , LOS 
NUE 

STOP 

END 

THIS  SUBROUTINE  CREATES  THE  INDEX  ARRAY 
FOR  AN  INVENTORY  DATA  ELEMENT. 

SUBROUTINE  SUM  ( INDEX , I . J ,K , L , NUM ) 

INTEGER  INDEX. I. J.K,L,NUM 

INDEX  =  I +  J*  1000000  + K-;100000+L*  1000  + NUM 

WRITE  (11.300)  INDEX 

FORMAT  (19) 

RETURN 

END 


Figure  D.4  WATFIV  Program  INV 

should  be  character,  vice  numeric,  arrays,  and  CMSIO  allows 
this  choice.  The  APL  functions  INVMATX  and  MATRIX  discussed 
below  assume  the  APL  character  arrays  are  words  of  ten  char¬ 
acters,  the  first  nine  characters  being  the  data  index,  and 
the  tenth  character  a  blank. 

3.  CREATING  THE  INVENTORY  AND  LOSS  ARRAYS 

Using  the  INVXX  arrays  created  above,  and  the  APL  func¬ 
tions  GETINV  in  Figure  D.6  and  INVMATX  in  Figure  D.7,  create 
the  arrays  IXX.  Note  that  GETINV  calls  INVMATX,  and  INVMATX 
uses  the  INVXX  arrays.  APL  workspace  size  limitations  may  be 
a  problem  due  to  the  large  amount  of  data,  and  it  may  be 


86 


$  JOB 


INTEGER  GRADE  ,  or' ,  YEAR ,  LOS  ,  IYR ,  LOSS 

THE  DO  LOOPS  CORRESPOND  TO  THE  NUMBER  OF 
RECORDS  FOR  YEARS  1977  THRU  1983. 

DO  100  1=1,2203 
DO  100  1=1,2231 
DO  100  1=1,2337 
DO  100  1=1,2351 
DO  100  1=1,2317 
DO  100  1=1,2324 
DO  100  1=1.2330 

READ  ( 10 , 2u , END=200 )  YEAR , OF , GRADE , LOS , INV , 
*L1,L2:l3  L4,L5,L6,L7,L8 
FO&MA'f  74X2  1 9l5  ) 

IYR  =  YEAR- 76 
LOSS  =  100000000 "IYR 

IF  (L.1.GT.0)  CALL  SUM  ( IN  ,  LOSS  ,  OF ,  GRADE  ,  LOS 

*  L 1 ) 

iF  [L2.GT.0)  CALL  SUM  (IN , LOSS , OF .GRADE , LOS 

*  L2  ] 

iF  [L3.GT.0)  CALL  SUM  ( IN , LOSS , OF , GRADE , LOS 

a1f (l4 . GT . 0  )  CALL  SUM  ( IN , LOSS , OF , GRADE , LOS 

,iF  / L5 . GT . 0 )  CALL  SUM  ( IN , LOSS , OF , GRADE , LOS 

,lF  (L6.GT.0)  CALL  SUM  ( IN , LOSS , OF , GRADE , LOS 
L6  ) 

,iF  (L7.GT.0)  CALL  SUM  ( IN , LOSS , OF , GRADE , LOS 

.IF  (L8.GT.0)  CALL  SUM  ( IN .LOSS , OF , GRADE , LOS 
'f ,  L8  ) 

CONTINUE 

STOP 

END 

THIS  SUBROUTINE  CREATES  THE  INDEX  ARRAY 
FOR  A  LOSS  DATA  ELEMENT. 

SUBROUTINE  SUM  ( INDEX , I , J ,K ,L , NUM) 

INTEGER  INDEX , I , J  ,  K , L , NUM 

INDEX  =  I  +  J  "  10  0  0  0  0  0  +  K  "  100000  +  L'"  1000  +  NUM 

WRITE  (11,300)  INDEX 

FORMAT  (19) 

RETURN 
END 


Figure  D.5  WATFIV  Program  LOSS 


necessary  to  create  one  or  two  arrays  at  a  time,  and  then 
copy  them  to  another  workspace. 

The  LXX  arrays  are  created  in  a  manner  similar  to  the 
above,  using  the  APL  functions  GETLOSS  in  Figure  D.8,  and 
MATRIX  in  Figure  D.9.  MATRIX  uses  the  loss  arrays  LOSSXX. 


Figure  D.6  APL  Function  GETINV 


V  Z+INVMATX  X ; A ;B ; C : D :E : F ; I :J 
p  CREATES  THE  INVENTORY  ARRAYS  FOR  THE  FISCAL 
p  YEARS  USING  THE  ARRAYS  OF  INDEXES  INVXX. 
p  INVXX  MUST  BE  A  CHARACTER  VECTOR  OF  9  DATA 
p  ENTRIES  FOLLOWED  BY  1  BLANK  FOR  EACH  LOOP . 
Z<-(40  31  10  )p0 
I«-pX 

«/-*•(  1  + 1  )*10 
LOOP:+(JsO)/OUT 
p  A-<-$(ltX) 

B-«-l  +  ($(2+X«-(l+X)  )  ) 

C<-l  +  (a>(ltX<-(2+X))  ) 

D<-1  +  ($(2+X-«-(l+X  5 )  5 

E-*-  SB  (3tX<-(2't-X)) 

ZIB-,D;C1+E 
X-*-(  4+X) 

«7«-«7-  1 
■+LOOP 

OUT: 'FINISHED  --  SHAPE  OF  MATRIX  IS  ' 


Figure  D.7  APL  Function  INVMATX 


To  create  the  aggregate  (e.g.,  aviation,  combat  support, 
ground  combat)  inventory  arrays,  use  the  APL  functions  GETAV 
in  Figure  D.10,  GETCS  in  Figure  D.ll,  GETGC  in  Figure  D.12, 
and  GETOF  in  Figure  D.13.  GETOF  is  called  by  the  other  three 
functions.  Each  calling  function  creates  on  aggregate  array, 


V  GETLOSS 

[1]  p  THIS  FUNCTION  CALLS  MATRIX  FOR  EACH  FISCAL 

[2]  p  YEAR.  LXX  IS  THE  LOSS  ARRAY  FOR  FISCAL  YEAR 

[3]  p  XX  BY  OF/LOS/GRADE. 

[4]  L77+MATRIX  L0SS77 

[5]  L7 8+MATRIX  LOSS 78 

[6]  L7 9+-MATRIX  LOSS 7  9 

[7]  L80+MATRIX  LOSS 80 

[8]  p  L81+MATRIX  LOSS8\ 

[9]  P  L82*-MATRIX  LOSS 8  2 

[10]  p  L83+MATRIX  LOSS 83 

V 


Figure  D.8  APL  Function  GETLOSS 


[1] 

[2] 

[3] 

[4] 
C5] 
C6] 

[7] 

[8] 

[9] 

[10] 
[11] 
[12] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[20] 
[21] 
[22] 


V  Z+MATRIX  X;A;B;C;D‘,E;F;I;J 
THIS  FUNCTION  CREATES  THE  LOSS  ARRAY  FOR 
THE  FISCAL  YEARS  USING  THE  ARRAY  OF  LOSS 
INDICES  LOSSXX .  IT  IS  CALLED  BY  GETLOSS. 
LOSSXX  MUST  BE  A  CHARACTER  VECTOR  WITH  9 
DATA  ENTRIES  FOLLOWED  BY  1  BLANK  FOR  EACH 
LOOP. 

Z*~  (40  31  10  )p0 
I+pX 

J<-(I  +  1)t10 
LOOP:+(Js 0 )/OUT 
( 1  fX  ) 

B*- l  +  (4(2+X-«-(14X)  )  ) 

C+- l  +  (a>(l+X<-(24-X)  )  } 

D-*-l  +  ($(2tX*-(l+Xj  )  ) 

E+vt  ( 1  +  X*-  ( 2  +X  )  ) 

F<-<&  (2+X<-vl+X;  ; 

ZIB :  D  :C]  +ZIB  ,D  ;C]  +F 
X*-  ( 3  4X  ) 

J  <-«/-  1 
+LOOP 

OUT: ' FINISHED  --  SHAPE  OF  MATRIX  IS' 

PZ 


Figure  D.9  APL  Function  MATRIX 


e.g.,  GETGC  creates  the  ground  combat  inventory  array  named 
GC.  Again,  workspace  size  limitations  may  be  a  problem,  and 
the  functions  can  be  altered  to  create  part  of  an  array  at  a 
time.  This  permits  the  operator  to  reduce  the  number  of  IXX 
arrays  present  in  the  workspace.  Note  that  the  arfays  are  a 


'central  inventory'  array,  taking  the  maximum  of  the  yearly 
loss  and  the  average  of  the  year's  inventory. 


7  GET  A  V 

THIS  CREATES  THE  INVENTORY  MATRIX  FOR 
THE  AVIATION  GROUP.  CALLS  GET OF . 
OCCUPATION  GROUP.  CALLS  GETOF . 

USE  LINES  6-7  FOR  ESTIMATION  INVENTORY 
MATRIX . 


AV*(n  31  10)p0 
AV-  GETOF  3  8 
R  USE  LINES  9-11 
n  MATRIX. 

V  AV*-  (  4  31  10)p0 
VAV-  GETOF  3  8 

V 


VALIDATION  INVENTORY 


Figure  D.10  APL  Function  GETAV 


V  GETCS 

a  THIS  CREATES  THE  INVENTORY  MATRIX  FOR  THE 
p  COMBAT  SUPPORT  GROUP.  CALLS  GETOF. 
p  USE  LINES  5-12  FOR  ESTIMATION  INVEN- 
p  TORY  MATRIX 
CS* ( 3  4  31  10)p0 
p  THIS  CREATES  THE  OF  ENGINEERS 
QS- \ • • • _ GETOF  7 

p  THIS' CREATES  THE  OF  OPERATIONAL  COMMUN- 
p  I CAT IONS 

2 • • • _ GETO F  13 

p  THIS’ CREATES  THE  OF  MOTOR  TRANSPORT 
CS- 3 ; ; ; -  - GETOF  2  0 

p  USE  LINES  15-22  FOR  VALIDATION  INVEN- 
p  TORY  MATRIX 
VCS* ( 3  4  31  10)p0 
a  THIS  CREATES  THE  OF  ENGINEERS 
VCS- 1 ; ; ; -  - GETOF  7 

R  THIS  CREATES  THE  OF  OPERATIONAL  COMMUN- 
r  I CAT IONS 
VCS-  2 • *  *  -  - GETOF  1 3 

r  THIS’ CREATES  THE  OF  MOTOR  TRANSPORT 
VCS-  3 : :  ;  -  GETOF  2  0 


Figure  D.ll  APL  Function  GETCS 


The  aggregate  loss  arrays  (e.g.,  AVL,  CSL,  GCL,  etc.) 
are  now  created  using  the  APL  functions  GCLOSS  in  Figure 


V  GETGC 

A  THIS  CREATES  THE  INVENTORY  MATRIX  FOR  THE 
p  GROUND  COMBAT  OCCUPATION  GROUP.  CALLS 
a  GETOF. 

A  USE  LINES  6-12  FOR  ESTIMATION  INVEN- 
A  TORY  MATRIX 
GC+( 3  4  31  10)p0 
A  THIS  CREATES  THE  OF  INFANTRY 
GC- 1 ; : ; - - GETOF  3 

A  THIS  CREATES  THE  OF  ARTILLERY 
GC- 2: ; i-- GETOF  5 

a  THIS  CREATES  THE  OF  TANKS  AND  AMPHIB 
GC- 3 ; ; ;  -  -  GETOF  1 0 

a  f/SF  LJiVES  15-21  FOF  VALIDATION  INVEN- 
a  TORY  MATRICES 
VGC+-  ( 3  3  31  10  )p0 
a  THIS  CREATES  THE  OF  INFANTRY 
VGC- 1 ; ; ;  -  -  GETOF  3 
a  THIS  CREATES  THE  OF  ARTILLERY 


VGC-  2 ; ; 


VGC-  3 : : 


--GETOF  5 


THIS  CREATES  THE  OF  TANKS  AND  AMPHIB 


--GETOF  10 


Figure  D.12  APL  Function  GETGC 


Z- GETOF  X;J;Y 

CREATES  THE  CENTRAL  INVENTORY  FOR  THE  YEAR 
DESIRED  USING  THE  ARRAYS  IXX. 


J+X  + 1 

^  USE  LINES 


Z<-(4 

31 

10  )c 

)0 

ZC1 ; 

K 

(Y*- 

1 77 

F  J; 

;3  +178 

r  J 

Z  [2  ; 

;  J-M 

(Y+- 

J78 

1  ,/; 

; 3  +17  9 

1  J 

Z  [  3  ; 

;  3-H 

U- 

179 

r  J ; 

; 3  +180 

i  J 

Z[4  ; 

:  3^< 

(Y<- 

ll  80 

Lc/ ; 

; 3  +18 1 

L  J 

USE  LINES  10-14  FOR  THE 
Z+(3  31  10)p0 
Z  Cl ;  ;]■«-((. Y <-(181  [J;  ;]+I82 


FOR  THE  ESTIMATION  MATRIX. 

CJ;  ;3+I78[t7;  ;]  U2)T  (L77CJ;  ;] 
[«/;  ;  3  +179  [</;  ;  3  2  )  i  (£78[</;;3 

[J;  ;3+I80[e7;  ;3  )*2  )T (L7  9 [J;  ;3 
[J;  ;3+I81[J;  ;3  )*2  )T  (L80[«7;  ;3 
4  FOR  THE  VALIDATION  MATRIX . 


Z[2;  ;]<- 


I82[J; ; 3  +183 


)  J  '  \  V  \  J.  u  4.  LW  f  t  J  ’  ^ 

Z  [3  ;  j3«-I83[J;  ;3r£83[«7;  ;3 


[«7 ; ;  3  )*2)T  (t81[Ji 
[J;  ;  3 )  +  2 )f  (£82  [J  j  ; . 


Figure  D.13  APL  Function  GETOF 

D.14,  AVLOSS  in  Figure  D.15,  and  CSLOSS  in  Figure  D.16 
These  functions  call  the  arrays  LXX  directly.  Again 
workspace  size  may  require  creative  changes  to  th 
functions . 


V  AVLOSS 

fl  THIS  CREATES  THE  LOSS  MATRIX  FOR 
p  AVIATION. 

R  USE  LINES  4-8  FOR  ESTIMATION  MATRIX. 
AVL+( 4  31  10)p0 
AVLl  1;  ;3<-L77[39;  ;] 

AVLL 2;  ;  3 -«-L7  8  [39  ;  ;3 
A7L[3;  ;3-«-L79[39;  ;] 

A7L[4;  ;  ]<-L8  0  [39  ;  ;3 

p  USE  LINES  10-13  FOR  VALIDATION  MATRIX 
VAVL+i  3  31  10)p0 
IM7LC1;  ;  ]  *-L8 1  [3  9  ;  ;] 

VAVLL 2;  ;]•*■£ 82C39;  ;3 
VAI7LE3;  ;  3-<-L83  [39  ;  ;] 

V 


V  CSLOSS 

p  FUJS  CREATES  THE  LOSS  MATRIX  FOR  COMBAT 

r  SUPPORT. 

r  USE  LLUES  4-16  FOE  ESTIMATION  MATRIX . 
CSL*- ( 3  4  31  10)p0 
CSL[1 ;  1 ;  ;3+L77[8;  ;] 

CSL  [  1 ;  2  ;  ;  3<-L7  8  [8  ;  ;  ] 

CSL [ 1 ; 3 ; ;3+L79[8; ;] 

CSS  [  1 ;  4  ;  {]«-£  80[8;  ;3 
CSL  [  2  ;  1  ;  ;]^-L77[14;  ;] 

CSL  [2  ;  2  ;  ;  ]«-£78C14;  ;] 

CSL [ 2 ; 3 ; ;]^L79[14; 

CSLl  2  ;  4  ;  ;  ]-(-L8  0  [14  ;  ;  ! 

CSL  [  3  ;  1 ;  ;3«-L77[21;  ;3 
CSL  [3  ;  2  ;  ;  3<-L78  [21 ;  ;  ] 

CSL  [ 3  ;  3  ;  ;  ] -e-L7 9  [  2 1  ;  ;3 
CSL  [3  ;  4  ;  ;  ]  *-L80  [21  ;  ;3 

fl  USE  LINES  18-27  FOR  VALIDATION  MATRIX. 
VCSL+(3  3  31  10)p0 
7CSL[l;l; ;3<-L81[8; ;] 

VCSL  [  1 ;  2  ;  ;3<-L82[8;  ;] 

FCSLC1 ;  3  ;  ; 3+L83  [8  ;  ;  3 
7CSLC2;1;  ;3<-L81[14;  ;3 
^CSL[2;  2;  ;3<-L82[14;  ;3 
VCSL [ 2 ; 3 ; ;3+L83[14; ;3 
l^CSL  [ 3  ;  1 ;  ;3«-L81  [21 ;  ;  3 
VCSL  [  3  ;  2  ;  ;  3  “-L82  [21 ;  ;3 
VCSL [3  ;  3  ;  ; 3+L83  [21 ;  ;  3 


VCSL [ 1 
VCSL [ 1 
VCSL [ 1 
VCSL  [  2 
VCSL  [  2 
VCSL [ 2 
VCSL  [  3 
VCSL [ 3 
VCSL [ 3 


V  GCLOSS 

a  THIS  CREATES  THE  LOSS  MATRIX  FOR  GROUND 

a  COMBAT. 

p  USE  LINES  4-16  FOR  ESTIMATION  MATRIX. 
GCL+ (3  4  31  10)p0 
GC£[1;1;  ;3-<-£77[4;  ;] 

GCL  [  1 ;  2  ;  ;3-*-£78[4;  ;] 

GCL  [  1 ;  3  ;  ;3-*-£79[4;  ;] 

GCL  [  1 ;  4  ;  ;3-*-£80[4;  ;] 

GC£[2;1;  ;  ]«-L77  [6  ;  s  ] 

GCL [2 ; 2 ; ;3^£78[6; ;] 

GCLi  2;3;  ;  ]  •*-£79  [6  ;  ;] 

CC£[2;4;  ;]-*-£80[6;  ;3 
GCL C 3  ;  1 ;  ;3-*-£77[ll;  ;3 
GCLI  3  ;  2  ;  ;  3 -*-£78  [11;  ;] 

GCL  [  3  ;  3  ;  ;  ] -*-£79  [  11 ;  ;3 
GCL  [  3  ;  4  ;  ;  ] -*-£80  [  1 1 ;  ;3 

a  USE  LINES  18-27  FOR  VALIDATION  MATRIX. 
VGCL+( 3  3  31  10)p0 


VGCLC 1 

fcc£  C 1 

VGCL [ 1 
FCC£[2 
FCC£  C  2 
VGCL 12 
VGCL [3 
VGCL  1 3 
VGCL [3 
V 


3-*-£81[4;  ;] 

3 -<-£8 2  [4  ;  ;  3 
3 -*-£83  [4 ;  ; 
3-*-£81  [6  ;  ; 
3-*-£82  [6  ;  ; 

3  -*-£8  3  [ 6  ;  ; . 

]  -<-£81  [11;  ;  3 
3  -«-£8  2  [11;  ;  3 
3 -*-£83  [11 ;  ;  3 


Figure  D.16  APL  Function  GCLOSS 

Use  the  APL  functions  GETNA  in  Figure  D.17,  GETNC  in 
Figure  D.18,  and  GETNG  in  Figure  D.19  to  create  the  desired 
grade  specific  inventory  arrays,  e.g.,  NA5  ,  which  is  the 
central  inventory  array  for  aviation  First  Lieutenants 
(grade  code  5).  See  Table  17  for  grade  codes.  These  func¬ 
tions  call  the  aggregate  inventory  arrays,  i.e.,  AV,  CS ,  and 


The  APL  functions  GETYA  in  Figure  D.20,  GETYC  in  Figure 
D.21,  and  GETYG  in  Figure  D.22  can  now  be  used  to  create  the 
desired  grade  specific  loss  arrays,  e.g.,  YA5 ,  which  is  the 
loss  array  for  aviation  First  Lieutenants.  These  functions 
use  the  aggregate  loss  arrays,  i.e.,  AVL,  CSL,  and  GCL. 

Use  the  APL  functions  GETANA  in  Figure  D.23,  GETANC  in 
Figure  D.24,  and  GETANG  in  Figure  D.25  to  create  the  arrays 


[1] 

[2] 

[3] 

[4] 

[5] 

C  6  3 

[7] 

[8] 

[9] 

[10] 
[11~ 
[12 

[13] 

[14] 

[15] 

[16] 
[17 
[18 
[19 
[20 
[21] 
[22 
[23 
[24 


V  GETNA 

a  THIS  CREATES  THE  GRADE  SPECIFIC  INVEN- 
q  TORY  FOR  AVIATION. 


a  USE  LINES  4-13  FOR  ESTIMATION  MATRICES. 
NA0+C3  1  2)5(1  4  31)pd7[;;l] 

Ml<-(3  1  2)15(1  4  31)pA7[;;2] 

NA  2«-(3  1  2  16}  C 1  4  31  )pAVt ;  ;  3] 

M3<-  3  1  2  1 6}  <  1  4  31)pA7[;;4; 

W4+  3  1  2)5(1  4  31)pA7[;;5] 

NA  5^(3  1  2)5(1  4  3lJpA7[;;6] 

NA  6<-(3  1  2)5(1  4  31  )pA7[ ;  ;  7] 

W7<-  3  1  2)5(1  4  31  )pA7[ ;  ;  8] 

AM8<-(3  1  2)5(1  4  31  )pA7[ ; ; 9] 

NA  9-6-  ( 3  1  2)5(1  4  31  jpAVt  ;  ;  10] 
a  USE  LINES  15-24  FOR  VALIDATION  MATRICES. 
VNA0+(3  1  2)5(1  3  3  1  )p7A7[ ;  ;  1] 

VNA  l<-(3  1  2)5(1  3  31  )p7A7[  ;  ;  21 
VNA 2^(3  1  2)5(1  3  3l)p7A7[;;3] 

7/VA3-e(3  1  2)5(1  3  3 1 )  p  7A  7  [  ;  ;  4  ] 

VNA 4^(3  1  2)5(1  3  31  )p7A7[;  ;  5] 

VNA5+13  1  2)5(1  3  31)p7A7[;;6] 

VNA6+C3  1  2)5(1  3  31  )p7A7[  ;  ;  7] 

VNA7+(3  1  2)5(1  3  3  1  )p7717[  ;  ;  8] 

VNA  8<-(3  1  2)5(1  3  31  )p7A7[  ;  ;  9] 

77VA9«-(3  1  2)5(1  3  31  )p7717[  ;  ;  10] 

V 


Figure  D.17  APL  Function  GETNA 

of  grade  specific  average  inventory.  For  example,  ANA5  is 
the  average  inventory  array  for  aviation  First  Lieutenants. 
The  above  three  functions  call  GETAAV  in  Figure  D.26,  GETACS 
in  Figure  D.27,  and  GETAGC  in  Figure  D.28,  respectively. 
Also,  this  second  group  of  three  functions  all  call  GETAOF 
in  Figure  D.29,  which  in  turn  uses  the  arrays  IXX . 

The  procedures  have  now  produced  central  inventory  and 
loss  arrays  for  the  estimation  and  validation  years,  and  the 
average  inventory  arrays  for  the  estimation  years.  All  of 
these  arrays  are  needed  to  calculate  attrition  rates  and 
figures  of  merit.  Note  there  is  no  requirement  to  create 
arrays  of  the  average  losses. 


V  GETNC 

[1]  A  THIS  CREATES  THE  INVENTORY  FOR  COMBAT 

[2]  a  SUPPORT. 

[3]  a  USE  LINES  4-13  FOR  ESTIMATION  MATRIX. 

[4]  NCQ*-  (.  3  1  2  )§CS[  ;  ;  ;1] 

[5]  NCl*-{  3  1  2  }§CS[  ;  ;  ;  2] 

[6]  NC2+- 1 3  1  2)S)CS[;;;3] 

[7]  WC3-e-(3  1  2  )S)CS[  ;  ;  ;4] 

[8]  ^4-f-(3  1  2  ISCSI;;;  5! 

[9]  NC5fr  (  3  1  2  )S)CS[;  ;  ;6] 

[10]  NC6+I3  1  2  )<s>CS [  ;  ;  ;  7 ] 

[11]  NC7+(3  1  2)SCS[;;;8] 

[12]  iVC8-e-(3  1  2  mcSLx ;  ;9] 

r  1 3  n  a ICQ'S-  ( 3  1  2  j r  •  •  •  1 0 1 

[14]  A  USE^  LINES  15- 24 ' FOR^ VALIDATION  MATRIX. 

[15]  YWC0*(3  1  2)SVCSC;;;1] 

[16]  VNCl<r\3  1  2  5 SPSS [  ;  ;  ;  2] 

[17]  VNC2+C3  1  2  WVCSZi  ;  ;3] 

[18]  VNC3  +  1 3  1  2 )SFCS[ ; ; ;4] 

[19]  VNCH+13  1  2  )ts>VCS[  ;  ;  {  5] 

[20]  VNC5+- 1 3  1  2  )SPCS[;  ;  ;6] 

[21]  VNC 6^(3  1  2)SFCS[;;;7] 

[22]  VNC 7-S-C3  1  2)SFCS[;;;8] 

[23]  FWC8*(3  1  2  )S)VCS[  ;  ;  ;  9] 

[24]  VNCQ+  ( 3  1  2 )SFCS[ ; ; ;10] 

V 


Figure  D.18  APL  Function  GETNC 


4.  ATTRITION  RATE  AND  FOM  CALCULATIONS 

The  remainder  of  the  procedures  assumes  certain  global 
APL  variables  are  defined  in  the  workspace: 

(1)  N  =  the  estimation  period  central  inventory  array, 
e.g. ,  NA5 , 

(2)  Y  =  the  estimation  period  loss  array,  e.g.,  YA5 , 

(3)  VN  =  the  validation  period  central  inventory 
array , e . g . ,  VNA5 , 

(4)  VY  =  the  validation  period  loss  array,  e.g.,  VYA5 , 

(5)  AN  =  the  estimation  period  average  inventory  array, 
e.g.,  ANA5 ,  and 

(6)  G  =  the  James-Stein  forced  shrinkage  rat^. 

All  functions  are  invoked  by  entering  the  function  name. 
Since  the  above  global  variables  and  the  global  output  vari¬ 
ables  from  functions  discussed  below  are  used  to  calculate 
attrition  rates  and  FOM's,  care  must  be  taken  to  ensure 
variable  values  are  not  changed  inadvertently. 


AD-A168  782 
UNCLASSIFIED 


LIMITED  TRANSLATION  SHRINKAGE  ESTIMATION  OF  LOSS  RATES  2/2 
IN  NARINE  CORPS  HANPOHER  NODELS(U)  NAVAL  POSTGRADUATE 
SCHOOL  HONTEREV  CA  J  R  ROBINSON  HAR  86 

F/G  3/9  NL 


V  GETNG 

fl  THIS  CREATES  THE  GRADE  SPECIFIC  INVEN 
a  TORY  FOR  GROUND  COMBAT. 
fl  USE  LINES  4-13  FOR  ESTIMATION  MATRIX. 
NG0+(3  1  2)<S)CCC;  ;  ;1] 

NG l<-(3  1  2}<S)GCl;  ;  ;2] 

NG2+- 1 3  1  2)5)CC[;;;3] 

NG3  +  13  1  2)i?GC[;;;4] 

WG4«-(3  1  2)S)CC[;  ;  ;5] 

NG5+C3  1  2  ;  ;6] 

NG 6+(3  1  2  )§GC C ; ; ;7] 

NG7+-  (  3  1  2  )S)CC[  ;  ;  ;  8] 

NG8+-  (  3  1  2)$GCC;;;9] 

NG9-*-  (  3  1  2  )SJCC[  ;  ;  ;  10] 
b  USE  LINES  15-24  FOR  VALIDATION  MATRIX 


VNG  0< 
VNG 1- 
VNG2* 


[3  1  2)«)7CC[ 
3  1  2)5?F£C[ 


VNG2*-(  3  1  2)S)7CC[ 
VNG3<rl3  1  2  ;($7CC [ 
7WG4<-(3  1  2)fS)FGC[ 
7/VC5-<-(3  1  2)$7CC[ 
W/C6-K3  1  2)S}FCC[ 
7WG7-*-(3  1  2)iS!7GC[ 
VNG8-t-(3  1  2)S>'/CC[ 
7M79<-(3  1  2  ;  [ 


Figure  D.19  APL  Function  GETNG 

Use  the  APL  function  ESTIM  in  Figure  D.30  to  calculate 
the  array  R  of  transformed  scale  attrition  rate  estimates. 
This  function  calls  BINPREP  in  Figure  D.31,  SUMSQ  in  Figure 
D.32,  and  MLE  in  Figure  D.33.  The  first  page  of  R  is  the 
array  of  original  scale  aggregate  rates,  each  entry  being 
equal.  The  second  page  is  the  array  of  transformed  scale 
aggregate  rates.  The  third  page  is  the  cell  MLE  rates.  The 
fourth  page  is  the  TSCA  estimates.  The  fifth  page  is  the 
James-Stein  estimates,  and  the  sixth  and  subsequent  pages 
(if  any)  are  the  limited  translation  James-Stein  rate 
estimates  for  different  values  of  d. 

The  user  now  has  several  options.  First,  the  relative 
savings  loss  for  different  values  of  d  can  be  calculated 
using  the  APL  function  RELS  in  Figure  D.34.  Enter  ESTIM  to 
select  the  d  values,  ensuring  the  shape  of  R  conforms  by 
adding  enough  limited  translation  pages  to  handle  all 


V  GETYA 

[1]  P  THIS  CREATES  THE  GRADE  SPECIFIC  LOSS 

[2]  P  ARRAYS  FOR  AVIATION. 

[3]  p  USE  LINES  4-13  FOR  ESTIMATION  MATRICES 

[4]  Y40<-(3  1  2  )$  ( 1  4  31  }pAVLl  ;  ;  1] 

[51  YA 1-K3  1  2)S)(1  4  31  WVX[;  ;2] 

[6]  Y42<-(3  1  2  5<S)  C 1  4  3l)pAVLCi  i3l 

[7]  YA3<-(3  1  21SK1  4  31  )pAVX[  ;  ;4] 

[81  YA4«-(3  1  2)S)(1  4  3lWPL[;;5] 

[9]  YA5<-(3  1  2)<9(1  4  3lWVL[;;6] 

[10]  YA6t-(3  1  2)S)(1  4  31  )pAVLl ;  ;  7] 

[11]  YA7<-(  3  1  2  )  5)  ( 1  4  31  ipAVLL  ;  ;  8] 

[12]  YA8«-(3  1  2  1  «S)  C 1  4  31  )pi47Z/[  ;  ;  9] 

[13]  YA9<-(3  1  2)19(1  4  31  )pAK£[; ;  10] 

[14]  p  USE  LINES  13-24  FOR  VALIDATION  MATRICES 

[15]  FYA0<-(3  1  2)19(1  3  31)pV71KL[;  ;1] 

[16]  FYAl-e(3  1  2)i9(l  3  31}pVAVLl;  ;2] 

I  [17]  VYA2+13  1  2)(9(1  3  3  1  )  p  VAVL  [  ;  ;  3  ] 

[18]  yYA3«-(3  1  2)19(1  3  3 1  )p  VAVL  [  ;  ;  4] 

[19]  VYAm-(  3  1  2)i9(l  3  31  }pVAVLl ;  ;  5] 

[20]  7YA5^(3  1  2)(9(1  3  3 1  )  p  VAVL  [  ;  ;  6  ] 

[21]  VYA6«-(3  1  2)i9(l  3  3 1  )p  VAVL  [  ;  ;  7] 

|  [22]  VYA7+C3  1  2  Sj(l  3  3 1  )  p  VA  VL  [  ;  ;  8  ] 

[23]  UY48«-(3  1  2)i9(l  3  3 1  )oVAVL  [  ;  ;  9] 

i  [24]  7Y^9^(3  1  2)i9(l  3  3 1  )p  VAVL  [  ;  ;  10  ] 

I  V 


Figure  D.20  APL  Function  GETYA 


selected  values  of  d.  Next,  enter  RELS  and  ensure  the  shape 
of  the  variable  RSL  (which  is  the  outputed  relative  savings 
loss)  conforms  to  the  number  of  d  values.  Then  run  RELS.  A 
graph  of  d  versus  RSL  allows  a  choice  of  the  best  d. 

The  second  option,  if  a  value  of  d  is  available,  is  to 
calculate  the  risks  of  the  various  estimation  methods  (in 
the  transformed  scale)  represented  in  the  R  array.  Use  ESTIM 
and  the  APL  function  RISKT  in  Figure  D.35  for  these  calcula¬ 
tions.  ESTIM  calculates  the  R  array  using  a  single  d  value, 
and  RISKT  uses  this  R  to  calculate  risks.  The  function  can 
only  handle  an  R  array  that  has  one  page  of  limited 
translation  rates.  Output  is  by  validation  year. 

Calculation  of  the  risk  in  the  original  scale  is  the 
third  option.  Use  the  function  RISKO  in  Figure  D.36  just  as 
RISKT  was  used.  This  function  calls  BINCONV  (Figure  D.37). 
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V  GETYC 

a  THIS  CREATES  THE  GRADE  SPECIFIC  LOSS 
a  ARRAYS  FOR  COMBAT  SUPPORT. 
a  USE  LINES  4-13  FOR  ESTIMATION  MATRICES 


YCO-H 

3 

1 

2 

GlCSLl 

9 

1] 

YC  1<- 

3 

1 

2 

§CSL  [ 

9 

2] 

YC2-H 

3 

1 

2 

GiCSLL 

9 

3] 

YC  3-e 

3 

1 

2 

§CSL  [ 

9 

4] 

YC4-H 

3 

1 

2 

GlCSLC 

• 

9 

5! 

YC5-H 

3 

1 

2 

SCSLL 

• 

9 

6] 

YC6-H 

3 

1 

2 

S )CSLl 

• 

9 

7] 

YC7-H 

3 

1 

2 

§CSLl 

9 

8] 

YC8<-< 

3 

1 

2 

S!CSL[ 

• 

9 

9] 

YC9«- 

3 

1 

2 

S)CSL[ 

9 

10] 

a  USE  LINES  15-24  FOR  VALIDATION  MATRICES 


VYC0+ 
VYC1  + 
VYC2  + 
VYC3  + 
FYC4<e 
VYCS  + 
VYC6+ 
VYC7*- 
VYCQ<r 
VYC  9<- 
V 


(312; 

§VCSL [ 

(3  1  2 

S)tY7SL[ 

(3  12 

StVCS  LL 

(3  12 

GlVCSLl 

3  1  2 

GiVCSLl 

(3  12 

SVCSLl 

(3  12 

GiVCSL  [ 

(312 

S*W7SL[ 

(3  12 

GiVCSLl 

(3  1  2] 

GlVCS  Lt 

; ;  l] 

; ; 2] 

; ;  3] 

;  ;4] 

; ;  5] 

; ; 6] 

; ; 7] 

; ; 8] 

; ;  9] 

; ;  10] 


Figure  D.21  APL  Function  GETYC 

Both  RISKT  and  RISKO  require  an  array  SS  to  be  defined 
globally.  SS  allows  selection  of  the  specific  cells  to  be 
included  in  the  risk  outputed  by  the  two  functions.  If  the 
global  risk  is  desired,  SS  should  be  an  array  of  all  ones. 
If  certain  cells  are  to  be  investigated,  e.g.,  cells  with 
inventory  less  than  six,  set  up  SS  so  that  each  cell  with  an 
inventory  less  than  six  has  an  entry  of  1  in  SS,  and  all 
other  cells  have  an  entry  of  zero.  SS  must  have  the  same 
shape  as  the  inventory  and  loss  arrays,  e.g.,  NA5  and  YA5 . 

The  fourth  option  is  to  calculate  the  original  scale 
attrition  rate  estimates  using  the  APL  function  BINCONV  in 
Figure  D.38.  Use  the  array  R  and  the  central  inventory  array 
for  the  aggregate  being  investigated. 


V  GETYG 

[1]  o  THIS  CREATES  THE  GRADE  SPECIFIC  LOSS 

[2]  o  ARRAYS  FOR  GROUND  COMBAT  COMBAT. 

[3]  o  USE  LINES  5-14  FOR  ESTIMATION  MATRICES 

[4]  YGO+  ( 3  1  2  )S)GCL[  ;  ;  ;  1] 

[51  YCl<-(3  i  2  )§GCL  [  ;  ;  ;  2  J 

[6]  YG2+(3  1  2  )(s)GCLC;  ;  ;3J 

[7]  YG3+ l 3  1  2 )§GCL [ ; ; ; 4] 

[8]  YC4W3  1  2  )<s>GCLl ;  ;  ;  5] 

[9]  YC5<-(3  1  2)$CCL[;;;6] 

[10]  YC6«-(3  1  2  )$}CCL[  ;  ;  ;7] 

[11]  YG7+-  (  3  1  2)$GCX[;;;8J 

[12]  YG8+ ( 3  1  2 }§GCL [ ; ;  ;  9] 

[13]  YGQ+-  (  3  1  2  )§GCL  [  ;  :  ;  10] 

[14]  o  USE  LINES  16-25  FOR  VALIDATION  MATRICES 

[15]  7YC0<-(3  1  2  )<siVGCLl ;  ;  ;  1] 

[16]  VYG1+C3  1  2  )<S)VGCLt ;  ;  ;  2] 

[17]  VYG 2«-(3  1  2  )§VGCLi  ;  ;  ;  3] 

[18]  VYG3+C3  1  2  WCCL[ ;  ;  ;4] 

[19]  KYC4W3  1  2  }<S)VGCLl ;  ;  ;  5] 

[20]  VYG  5f(3  1  2  )<s>VGCLt ;  ;  ;6] 

[21]  VYG&+- ( 3  1  2)$7CCL[js;7] 

[22]  YYC7«-(3  1  2  )«)7<;CL[;  ;  ;8] 

[23]  UYC8<-(3  1  2  WCCL[  ;  ;  ;  9] 

[24]  7YC9-*-(3  1  2  )§VGCL  [  ;  ;  ;  1 0 ] 


V  GET AN A ; Z 

[1]  n  THIS  CREATES  THE  AVERAGE  INVENTORY  FOR 

[2]  o  AVIATION.  CALLS  GET AAV . 

[3]  fl  i/SE  LliVES  4-14  EOE  ESTIMATION  ARRAYS 

[4]  z<-cez7mi7 

[5]  4/V/?0<-(3  1  2  )§  ( 1  4  31)pZ[;;l] 

[6]  AM1<-  3  1  2  5«S}  C 1  4  31)pZ[;;2] 

[71  AM2<-  3  1  2)S)(1  4  31  )pZ [ ;  ;  3] 

[8]  ANA 3<-(3  1  2)SH1  4  31)pZ[;;4] 

[9]  M^4<-  3  1  2)<5(1  4  31  )pZ[ ; ;  5] 

[10]  3  1  2)<5(1  4  31)pZ[;;6] 

[11]  ANA6+C3  1  2)<5(1  4  31)pZ[;;7] 

[12]  ANA7<r(3  1  2)t5(l  4  31)pZ[;;8] 

[13]  ANA8<r(3  1  2)<5(1  4  31)pZ[;;9] 

[14]  ANA9+(3  1  2)S)(1  4  31  )pZ  [  : ; 10] 

[15]  o  USE  LINES  16-26  FOR  VALIDATION  ARRAYS 

[16]  Z +GETAAV 

[17]  7M40<-(3  1  2)$(1  3  3  1  )pZ  [  ;  ;  1] 

[18]  VANA1  +  C3  1  2)15(1  3  31)pZ[;;2] 

[19]  VAN  A  2+ (3  1  2)  (5(1  3  31)pZ[;;3] 

[20]  7Mi?3<-  3  1  2)S)(1  3  31)pZ[;;4] 

[21]  VANAH+13  1  2  mi  3  3 1  )pZ  [  ;  ;  5  ] 

[22]  VANA5+C3  1  2)15(1  3  31)pZ[;;6] 

[23]  VANA6+13  1  2)15(1  3  31)pZ[;;7] 

[24]  VAN  A  7<-(3  1  2)15(1  3  31)pZ[;;8] 

[25]  VANA%+\3  1  2)«5(1  3  31)pZ[;;9] 

[26]  VANA<3+{3  1  2)(5(1  3  31)pZ[;;10] 
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V  GETANC 

A  THIS  CREATES  THE  AVERAGE  INVENTORY  FOR 
A  COMBAT  SUPPORT.  CALLS  GET ACS . 

A  USE  LINES  4-13  FOR  ESTIMATION  ARRAYS. 
Z+GETACS 

ANC  0«-(3  1  2)$Z[;;;1] 

ANC1+13  1  2  )S)ZC  ;  ;  ;  2! 

ANC 2<r(3  1  2  )<9Z[  ;  ;  ;  3] 

ANC3<r(3  1  2)S3ZC;;;4] 

ANC1**- 1 3  1  2  )5)Z[  ;  ;  ;  5] 

ANC  5<-(3  1  2  J§Z  [  ; ; ; 6] 

ANC6*r(3  1  2  5 S}Z  C  ;  ;  ;  7 1 

ANC7*-(3  1  2  )S}Z[  ;  ;  ;  8] 

ANC  8<-(3  1  2)S)ZC;;;9] 

ANCQ*-  (  3  1  2  )S!Z[  ;  ;  ;  10] 


ANC2*r(3  l  2)i?ZC 
ANC3*r(3  1  2)S3Z[ 
AM74<-(3  1  2)$)Z[ 


AiVC5< 
ANC  6« 
A/VC7< 
AWC8< 
4WC9< 


,3  1  2)S)ZC 
,3  1  2  5 S}Z  C 
,3  1  2)S|Z[ 
,3  1  2)S)Z[ 
!  3  1  2)<9Z[ 


a  (/SJ?  LIWES  16-26  FOR  VALIDATION  ARRAYS. 


u  '  v/ui.  n  v. 

I/71/VCCN-< 

3  1  2 

«?zc 

WlNCl-H 

3  12 

$zc 

VANC2*- 

3  12 

S)Z  [ 

VANC3*- 

3  12 

fs?z  L 

VANC^*-{ 

3  12 

fs}Z  [ 

VANC  5*- 

3  12 

(9Z[ 

VANC6*- 

3  12 

$Z[ 

VANC7+( 

3  12; 

JsjZ  [ 

VANC8 *- 

3  12 

$Z[ 

VANCQ+ 

’312] 

$Z[ 

V  Z+GETACS; Z 

p  THIS  CREATES  THE  AVERAGE  INVENTORY 
p  MATRIX  FOR  COMBAT  SUPPORT. 
p  CALLS  GETAOF. 

R  USE  LINES  5-11  FOR  ESTIMATION  ARRAYS. 
Z+(3  4  31  10)p0 

-  THIS  GETS  THE  OF  ENGINEERS 
Z- 1 ; ; GETAOF  7 

-  THIS  GETS  THE  OF  OPERATIONAL  COMM 

2  _ 2 •  •  • _ GETAOF  1 3 

-  THIS' GETS  THE  OF  MOTOR  TRANSPORT 
2- 3 • • •« . GETAOF  20 

p  use' Lines  13-19  for  validation  arrays 

p  Z«-(3  3  31  10  )p0 

-  THIS  GETS  THE  OF  ENGINEERS 


p  Z-  1 ; 

-  THIS 
p  Z-  2; 

-  THIS 
p  Z-  3  : 


'.--GETAOF  7 

GETS  THE  OF  OPERATIONAL  COMM 
. _ GETAOF  1 3 

GETS  THE  OF  MOTOR  TRANSPORT 
'.--GETAOF  2  0 


Figure  D.27  APL  Function  GETACS 


[1] 
C2: 
[3! 

[4] 

[5] 

[6] 
c  7 : 
[8] 
[91 
[10] 
[11] 
[12] 

[13] 

[14] 


-  Z-GETAOF  X;J;Y;Z 

a  CREATE  THE  AVERAGE  INVENTORY  FOR  THE 
a  YEARS  DESIRED  FOR  A  SPECIFIC  OF. 

JfrX  +  1 

a  USE  LINES  X-XX  FOR  ESTIMATION  ARRAY . 
Z+- ( 4  31  10)p0 

Z[l;  ;]■«-(  J7  7[«7;;]  +J78  LJ  i  ;  ]  )*2 
Z [2  5  ;]<-(J78[J;  ;]+J79[J;  ;]  5*2 
Z[3  ;  ;  ]<-(J79  [J ;  ;]+I80[J;  ;]  )f2 
Z[4; ;]+(I80[J; ;]+I81[e7; :]  )*2 
a  C/SF  LINES  XX- XX  FOR  VALIDATION  ARRAY. 
A  Z<-(3  31  10)p0 

a  Z  [  1 ;  ;]-«-(J81[c7;  ; ] +J82 [J ;  ;]  )»2 
a  Z[2;  ;]*(J82[J;  :]+I83[J;  ;]  )*2 
a  Z  [3  ;  ;  ]  <-T 8 3  [J ;  ;  J 
V 


Figure  D.29  APL  Function  GETAOF 
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[1] 

[2] 

[3] 

[41 

[5] 

[6] 

[7] 

[8] 

[9] 

[10] 
[11] 
[12] 
[13] 
[  1 4  J 

[15] 

[16] 

[17] 

[18] 

[19] 

[20] 
[21] 
[22] 
[23] 
[24; 

[25] 

[26] 

[27] 

[28] 

[29] 

[30] 

[31] 

[32] 

[33] 

[34] 

[35] 

[36] 

[37] 

[38] 


1  P<rESTIM;I;S;U;Cl;ClK;M;ZB;  ZBBA',  AGO  \AGT\D\P 
CALCULATES  THE  AGGREGATE  (.ORIGINAL). 
AGGREGATE  (TRANSFORMED) ,  MLE  (ORIGINAL), 
TSCA  (TRANSFORMED) .  JAMES-  STEIN ,  AND 
LIMITED  TRANSLATION  ESTIMATERS  FOR  Z  AS 
SCREENED  BY  D.  THE  SCREENING  MATRIX  FROM 
THE  ESTIMATION  YEARS  AVERAGE  INVENTORY. 

EACH  FACE  OF  R  IS  A  SEPARATE  ESTIMATE  IN 
THE  ESTIMATOR  ORDER  LISTED  ABOVE.  CHANGE 
LINES  12,15,35,36,37  TO  CHANGE  DEE'S  USED. 


N 


1  +  2  xP 


Z+Y  BINPREP 
DEE*  0 
J-e  6 

D+(+/  3  12  <$AN* 0)*0 
R+ (6 , (1+pZ) )p0 

s+sumsq  z 

P+(+/ .Yx(pY)pD)*+/ ,N*(pN)pD 
AGO+( (0 . 5  +  (+/A/*l+pJV)  )*0 . 5 )x  lo' 

AGT+1+/ .ZB)++/ ,D 
ZBBA+ (pZB )pZBB 
o on  -  SHJ+G 
(  2  +  1 +  HLC  ) x  \  0<C 

fl  SHJ  IS  THE'JS  SHRINKAGE  FACTOR. 

OOTl-  SHJ<r(K-  3  )*  (2-  K-  K*M)xt/S 
U+ I (ZB- ZBBA)*SHJ*0 .5 

RZ i:  ‘  '  ~  ~ 

RZ  2:._  .  „  . 

/?  [  3  ;  ;  ]+-Y  MLE 
]+Z  B 

]<-DxZBBA+C*ZB-  ZBBA 

IS  ARRAY  OF  LTJS  SHRINKAGE  FACTORS. 
LL  :  Cl<-0(  1- SHL+(llDEEiU)xSHJ 
H [I ;  ;  ]-s-DxZBSA+ClxZB-  ZBBA 
o  USE  THE  FOLLOWING  LINES  IF  MORE  THAN  ONE 
o  VALUE  OF  DEE  IS  USED. 
o  I+I+l 
R  DEE+-DEE  +  0.2 
R  +LL*\(DEE< 1.6) 

V 


RZ  4 
fl[5 

R  StfL 


)<-(l±pR)pAGO 
]+"(l  +  p/w  p  AGT 


N 


Figure  D.30  APL  Function  ESTIM 


[1] 

[2] 

[3] 

[4] 

[5] 


7  Z+Y  BINPREP  N ; P 

a  PREPS  THE  FREEMAN- TUKEY  VERSION  OF  THE 
o  ARC  SIN  TRANS  FOR  BINOMIAL  DATA 
a  Y  I S _LOSSES ;  N  IS  INVENTORY 
Z  +  ~  1  o  l+2xy+AT  +  l 

Z<-0 . 5X  (  (o  .  5+A/)*0 . 5  )xZ+  lo“l+2x  (Y  +  l  )+/V  + 1 
V 


Figure  D.31  APL  Function  BINPREP 


V  X+SUMSQ  Z ;SSE 

p  CALCULATES  THE  SSE  AND  SSB  FOR  Z.  ALSO 
p  CALCULATES  THE  RLE  ( ZB )  AND  GRAND  MEAN 
p  OR  AGGREGATE  (Z SB),  SO:TE  DIRECTLY  FROM 
p  TRANSFORMED  DATA. 

K++/,D 

ZB+D*  ( +  /Z  ) tM*-1  +  p Z 
Z BB+(+/ ,ZB)*K 
X++/, (Z- (pZ)pZBS)*2 
X+(X-SSE),SSE«-  +  /,  (Z-  (pZ  )pZS  )*2 
7 


Figure  D.32  APL  Function  SUMSQ 


V  Z«-F  MEE 

p  CTlECEMrES  27/E  MLE  IW  TEE  ORIGINAL  SCALE 
o  AND  TRANSFORMS  IT  INTO  ARCSIN  SPACE. 
D<-(+/  3  12  0)*0 

Ml«-<  (p+/7)pE)x  (+/y  )*  +  /# 

Z<-C*  ((0.  5  +  (+//V)t1  +  p/V)*0 . 5  )x  1 o- 1 +2  *Ml 


Figure  D.33  APL  Function  MLE 


Cl] 

[2] 

■  3- 
[4] 

[5. 
[6! 
C7; 
[81 

[9] 

[10] 
[11] 
[12] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[20] 


VRELSIU] V 

V  RELS ; RM ; RJ ; RL ;  I ;  Z 

o  COMPUTES  THE  RELATIVE  SAVINGS  LOSS  ( RSL ). 
a  IN  THE  TRANSFORMED  SCALE.  VY  IS  THE  ARRAY 
-  OF  ACTUAL  LOSSES.  VN  IS  THE  ARRAY  OF 
q  CENTRAL  INVENTORY.  R  IS  FROM  ESTIM.  VN . 
q  VY,R  MUST  BE  IN  THE  WORKSPACE.  CHANGE  LINE 
q  12  FOR  DIFFERENT  NUMBER  OF  DEE'S.  EACH 
q  FACE  OF  R  IS  THE  ESTIMATE  FOR  A  DIFFERENT 
q  DEE .  CALLS  BINPREP . 


q 


1+ 1 
Z+VY 
RSL 
RT 
RJ- 
LL  :RQ 


BINPREP  VN 
Z  )  )p0 
PZ'  '  ' 

,  PZ 
(1  +  5 


1  o  1  N  cl 

■+  (  9  ,  (1- 
+  /((l+f 
+  /(  (1  +  f 

Q+RZU- 


RL+  +  / ( ^ 1  +  d  Z 
RSLII\  ]<-(££ 
I+I  + 1 


] 


pZ  )pi?  [4 
p  Z  )  pi?  [  5  : 


Iix/(l+pZ)))p, (Z- (pZ)pRQ)*2 
RJ )tRT- RJ 


+LL*\  ((I  +  6)51+p/?) 

RSL 

V 


Figure  D.34  APL  Function  RELS 


7  RISKT ;  I ;  D ;  RM ;  RJ RL ;  RAO ;  RAT  \  RT ;  Z 
COMPUTES  THE  SCALED  RISK  IN  THE  TRANSFORMED 
SPACE.  VY  AND  VN  ARE  THE  VALIDATION  YEAR 
LOSSES  AND  INVENTORY.  VN ,  VY ,  R,  AND  AG 
MUST  BE  IN  THE  WORKSPACE. 


Z+VY  BINPREP  VN 
1+1 

RAO+RAT+RM+RT+RJ+RL+S  p  0 
ZK(+/  3  12  5MjV*0)*0 
LLiRAOUI  +  l*/  SSII-.  ;lx(Z 
RATin  +  l+Z  .SSII;  :  UlZLI 


:  ;]x(ZCJ:  ;]-/?[! 

RAFLU+1+/.SSII; ;  UlZLI;  ;  J  -  i?  [  2  ;  :1  )*2)  ++/.fi' 
i?M[  J]  •<-  ( +  /  ,  SSII  ;;Jx(Z[I;;j-/?[3;  ;]  )*2  )*+/,£> 
RTII1  +  1  +  /  ,SS[I;  ;  WZ[I;  ;]-i?[4;  ;]  *2  H+/,D 
i?e7[I]<-(  +  /,SS[I;  ;]x?ZCJ;  ;]-i?C5;  ;]  )*2U+/  ,D 
RLCI1  +  1+/ ,SSCI‘,  ;]*(ZCI;  ;]-R[6;  ;]  )*2  U  +  /  ,D 
I+I  +  l 

+LL* i (IS 3  ) 

•  tfIStf  I7V  TRANSFORMED  SCALE  USING  THE  ESTIM-  ' 
' ATORS  LISTED  BELOW.' 

'  81  82  83’ 


: ] )*  2 )r  +  /  ,D 
)*2 ) r+/,D 
*2 )r+/ ,D 
*2 )t  +  / ,D 
*2 )t  +  /  ,D 
*2  U  +  /  ,D 


81 

AGG  ORIG  RISK 

'  .  (jMO 

AGG  TRANS  RISK 

' , (vRAT 

MLE  RISK 

'  .  ( $RM ) 

TSCA  RISK 

'  ,  ($RT) 

JS  RISK 

<  ,  (® RJ) 

LTJS  RISK 

' ,  Url) 

VRISKO LQJV 

V  RISKOiDxRl;R2:R3',K:VxB;ARiS;R^\R5iRB;A;RR 
n  COMPUTES  THE  RISK  IN  THE  ORIGINAL  SPACE. 

o  yy  awd  iw  A/?ff  rw£:  wiLiMrraw  YffAff  losses 

ft  AWD  INVENTORY.  VN,VY,R.N  MUST  BE  IN  THE 
o  WORKSPACE.  CALLS  BINCONV . 

K+ 1 

£>-«-v/AW>0 

y<-  ( ( p  v )  p  d  )  x  v+binconv 
(pyy)pz? 
s«-(y*o  )A(y*i ) 

Rl+R2+R3<-Rm-RS<-R8<-  ( 1  +  pW  )pO 
£M:A<-(py)pAff[ff;  :3 

/?#<-(  (py)p£>)xsx(  (Py)Pw)x  (sx(^-  y)*2  )*/x(i-  y) 
R1CKJ++/ ,SSCKs  tJxRRCl: ;3 
ff2  [ff] <-  +  /  , SSCff ;  ;]x/?/?[2;  ;3 
KK«-+/  ,£> 

,S[3  ;  :  3 

R3  [K]<-+/ ,SSLK;  ;]*RRC3;  ;]xKKiKKl 
RHCKJ++/ ,SSCK; ;JxRR[n: ;3 
R5CK1++/ ,SSCK; ;]*RR[5; ;3 
ff6[ff3^+/,SS[ff; ;3*ffff[6; ;3 
K+K+l 

+LM*x (KSl+pVN) 

'RISK  IN  ORIGINAL  SCALE  USING  THE  ESTIMATORS 
'LISTED  BELOW.' 

'  81  82  83' 

'AGG  ORIG  RISK  =  '  ,  UR  1) 

'ACG  TRANS  RISK  =  ',UR 2) 

' MLE  RISK  5  '  , 3  ) 

'TSCA  RISK  =  ' , 

'JS  RISK  =  ',($/? 5  ) 

'LTJS  RISK  e  '  .  UR 6) 


9 

:®/?i ) 

f 

5/?2  ) 

9 

®/?3  ) 

9 

9 

Ur  5) 

• 

I$ff6  ) 

Figure  D.36  APL  Function  RISKO 


[•] 


V  B+BINCONV;VO-,V1-,N1;D 

Cl]  fi  INVERTS  ARC  SIN  TRANSFORMATION.  R  IS  THE 

[2]  R  AMM  TO  BE  INVERTED.  N  IS  THE  CENTRAL 

[3]  R  INVENTORY  FOR  THE  ESTIMATION  YEARS.  AN  IS 

[4]  p  THE  AVERAGE  INVENTORY  FOR  THE  ESTIMATION 

[5]  p  YEARS.  R,  N,  AND  AN  ARE  GLOBAL  VARIABLES, 

[6]  p  AND  MUST  BE  IN  THE  WORKSPACE. 

[7]  D<-v/AN>  0 

[8]  St-0. 5*l  +  lo/?f  <0. 5+(Nl*-(pR)p  +/N  )-rl  +  pAO*0. 5 

C9]  70t-H<  (- ot2  )x  (/Vl+0 . 5  )*0 . 5 

[10]  Vl<-R>  (ot2  )x(W1+0. 5  )  *  0  .  5 

[11]  B<-(  (pS)pS)xyi+Hxl/0^^1 

V 


Figure  D.37  APL  Function  BINCONV 
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APPENDIX  E 

RELATIVE  SAVINGS  LOSS 


1.  THEORETICAL  IMPLICATIONS  OF  RELATIVE  SAVINGS  LOSS 

The  value  of  d  can  be  empirically  chosen  using  a  concept 
developed  by  Efron  and  Morris  [Ref.  5]  called  the  relative 
savings  loss  (RSL).  The  RSL  quantifies  how  well  limited 
translation  does  versus  Jaines-Stein  estimation.  If 

(1)  R^.  =  risk  of  the  TSCA, 

(2)  Rj  =  risk  of  the  James-Stein  estimator,  and 


(3) 

then, 


RT  T  =  risk  of  the  limited  translation  James-Stein 
estimator , 


RSL  =  (Rlj-Rj)/ (Rt-Rj) .  (E.l) 

RSL  is  the  proportional  increase  in  global  loss  if  limited 
translation  estimation  is  used  instead  of  James-Stein 
estimation. 

Efron  and  Morris  [Ref.  5]  further  state  and  prove  a 
theorem : 

RSL  =  2[(d2+l)(l-cD(d))-d(p(d)]  (E  .  2  ) 

where  O  and  <p  are  the  standard  normal  c.d.f.  and  density 
function,  respectively.  By  this  theorem,  RSL  is  a  function 
of  d  only.  Figure  E.l  graphs  values  of  the  RSL  against  d. 

2.  APPLICATION  OF  RELATIVE  SAVINGS  LOSS 

Since  RLj  is  a  function  of  d,  a  vector  of  values  of  R^j 
can  be  calculated  for  d>0.  This  can  yield  a  vector  of 
values  of  the  RSL.  These  values  can  be  graphed  against  the 
values  of  d  which  generated  the  the  RSL  vector.  Efron  and 


01  80  90  to  20 

SS01  SONIAVS  3All\n3d 


Figure  E.l  Theoretical  Graph  of  RSL  versus 
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Morris  [Refs.  5,6]  state  that,  given  the  normality  with 
common  variance  assumption,  the  relationship  of  the  global 
risks  of  maximum  likelihood,  James-Stein,  and  limited 
translation  will  be 

RJ  <  RLJ  <  RT-  (E-3) 

Tables  4  thru  12  clearly  show  this  is  not  the  case  when  we 
identify  the  validation  FOM's  as  the  components  of  equation 
E.3. 

Figures  E.2  thru  E.7  are  the  graphs  of  calculated  RSL 
versus  d  for  the  study  cases.  The  encouraging  results  are 
misleading.  Tables  4,  5,  and  6  display  risk  values  that 
yield  negative  quantities  in  the  numerator  and  denominator 
of  equation  E.2.  Although  the  graphs  are  visually  correct, 
the  underlying  computations  are  at  variance  with  theory,  and 
are  simply  generating  offsetting  errors.  Therefore  RSL 
should  not  be  used  to  choose  d  for  the  present  aggregations. 


APPENDIX  F 

ROBUST  PARAMETRIC  EMPIRICAL  BAYES  ANALYSIS 


Robust  parametric  empirical  Bayes  (RPEB)  estimation  is  a 
version  of  a  generalized  linear  model  (see  McCullagh 
[Ref.  7])  combined  with  robust  Bayes.  It  attempts  to  fit  a 
probabilistic  model  that  assumes  the  cell  value  p^  for  the 
i^  cell  is  a  random  variable  that  can  be  described  with  a 
probability  distribution,  e.g.,  the  logistic  or  the  Poisson. 

Gaver  and  0 ’ Muircheartaigh  [Ref.  8]  developed  this  tech¬ 
nique  using  the  Poisson  to  model  collections  of  similar 
objects  that  independently  generate  events  in  accordance 
with  Poisson  processes.  First,  the  superpopulation  parame¬ 
ters  are  estimated,  using  point  and  interval  estimates. 
Next,  the  superpopulation  parameter  estimates  are  used  to 
compute  point  and  interval  estimates  of  the  individual  rate 
parameters.  The  result  is  a  Bayes  estimate  of  p^  ,  with 
limited  shrinkage  if  the  inventory  is  small. 

Gaver  has  suggested,  in  the  spirit  of  reference  8,  using 
a  logistic  model  for  p^  with  explanatory  variables  for  OF, 
LOS,  and  grade,  and  an  extra  binomial  variance  term  to 
partly  explain  differences  across  cells.  First,  using  numer¬ 
ical  integration,  the  maximum  likelihood  function  L(^,t)  is 
used  to  find  the  values  of  ^  and  t  that  maximize  the  likeli¬ 
hood  function.  Second,  the  expectation  of  p^  given  the  loss 
y L  is  calculated.  This  is  the  posterior  mean  of  p^,  and  is 
the  Bayes  estimator  of  p^,  with  limited  shrinkage, 
especially  if  the  inventory  n^  is  small. 

Therefore,  let 

Pi  -  em/ ( l+em)  =  e^/(l*eT)  (F.l) 


120 


where 


m  =  xi/3+5i 
q  =  M+T£i(n) 
i  =  1, . .  .  ,1 
5^  ~  rStudent  t 

=  losses  from  the  it^1  cell 

f'  K 

n^  =  inventory  in  the  icn  cell. 

Then  the  maximum  likelihood  function  is 

L(m,t)  =  n/{eM  +  rX/[(1  +  e^>TX)(1+(x/T)2/n)(n+l)/2])dx  (F.2) 

where  i  =  l,...,n,  and  the  integration  is  performed  over  the 
real  line.  This  likelihood  function  is  used  to  find  ^  and  r 
that  maximize  L(^,r)- 

Next,  if  we  assume  y ^ ,  and  let  y  =  y^,  then  the  expected 
value  of  p^  is 

E[Pi]  =  /pi(x)piy(l-pi)r  (F.3) 

c(n)/[h(x/T)J/n]^n+1^2dx, 

where  r  =  n  — y^,  c(n)  is  the  appropriate  normalizing 
constant,  and  the  integration  is  performed  over  the  real 
line.  This  expected  value  is  the  empirical  Bayes  estimator 
of  pj_. 

See  Deely  and  Lindley  [Ref.  9]  for  more  on  empirical 
Bayes . 
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